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A. Einleitung 
Die Fermische Idee und die Schwierigkeiten in ihrer Durchfihrung. 


Man kennt zur Zeit folgende Elementarteilchen: 


Bo und sein Antiteilchen 3*° 


Brae sy. 5s “ 70+ Hyperonen 

ue ae FS, A° 

BN? i et N-0 Nukleonen 

C30 eS . K~° Schwere Mesonen 

x*%- (x* ist Antiteilchen von x-; x° ist sein eigenes Antiteilchen) 
ut (u* ist Antiteilchen von p-) 

e*- —_ (e* ist Antiteilchen von e-) 

vv (v und ¥ sind verschieden, wenn Lepton-Erhaltung gilt) 
Y (ist sein eigenes Antiteilchen) 


Gelegentlich wird die Existenz weiterer Elementarteilchen diskutiert, entdeckt 


wurde bisher keines. 
Die Hauptsache ist, die Eigenschaften der bekannten Teilchen naher zu studieren 
und eventuell neue Teilchen zu finden. Diese Aufgabe sto8t auf Schwierigkeiten: 
— Erstens gibt es noch keine Theorie, die es gestatten wirde, die Existenz und 


die Eigenschaften gerade dieser Teilchen zu verstehen, 
— zweitens sind die meisten von ihnen instabil und man muB sie entweder in der 


kosmischen Strahlung beobachten oder sie kiinstlich erzeugen. In jedem Falle 
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entstehen sie bei ZusammenstéBen energiereicher Teilchen und dabei werden 
im allgemeinen gleich mehrere zugleich erzeugt. Das bedeutet, daf man es 
mit einem Mehrkérperproblem zu tun hat. 


Die Behandlung dieser Aufgabe ist aus zwei Griinden wiinschenswert : 


1. Man méchte die Vielfacherzeugung von Elementarteilchen in Hochenergie- 
prozessen um ihrer selbst willen verstehen. Dabei werden die bekannten Teil- 
chen mit ihren Eigenschaften als gegeben angesehen. 

2. Man méchte die Eigenschaften der einzelnen Teilchen und ihre Wechsel- 
wirkungen miteinander unter giinstigen Versuchsbedingungen beobachten. Da- _ 
zu braucht man Beschleuniger (Linearbeschleuniger, Cyclotron, Synchro- | 
cyclotron, Synchrotron), die energiereiche Primarteilchen auf — wenn moglich 
— Protonen schieBen. Dann finden die oben erwahnten Vielfacherzeugungs- 
prozesse statt und alle Sorten von Teilchen verlassen das Target (d.h. das 
beschossene Stiick Materie). Wenn die Theorie der Vielfacherzeugung zu be- 
rechnen gestattet, wieviele Teilchen jeder Sorte entstehen und welche Energie- 
und Winkelverteilung sie im Laborsystem haben, so gibt sie dem Experimental- 
physiker ein Mittel in die Hand, optimale Versuchsbedingungen fiir weitere 
Untersuchung zu schaffen. 


Selbstverstandlich schlieBt die Erledigung des ersten Punktes die des zweiten ein. 
Leider ist man noch weit entfernt davon, die Vielfacherzeugung restlos zu ver- 
stehen. Trotzdem wei man im Gebiet der Primarenergien zwischen 1 GeV und 
100 GeV soviel, da’ man wenigstens teilweise die beiden Anliegen erfiillen kann. 
Ein Mittel hierzu ist die statistische Theorie der starken Wechselwirkungen. 
Diese Theorie wurde erstmals von E. Fermi [J] vorgeschlagen, jedoch nach 
einigen (qualitativen) Anfangserfolgen als unzuverlassig beiseite gelegt. Erst 
nach einigen Jahren wurde sie wieder aufgegriffen und von mehreren Autoren so 
verbessert, daB sie heute verschiedene Fragen quantitativ beantworten kann. Die 
statistische Theorie beschreibt Produktion (und Vernichtung) derjenigen Ele- 
mentarteilchen, die starke Wechselwirkung untereinander haben. In der obigen 
Liste sind das alle von © bis x, wahrend die iibrigen Teilchen sowohl an die 
oberen als auch untereinander schwach oder elektromagnetisch gekoppelt sind. 
Wir werden uns auf die stark wechselwirkenden Teilchen beschranken. Ferner 
sollen hier nur Methoden besprochen werden, die im Gebiet von etwa 1—100 GeV 
(Lab) benutzt werden kénnen. Fiir extrem hohe Energien (> 103 GeV) verweisen 
wir auf den kiirzlich erschienenen zusammenfassenden Bericht [2]. 

Unsere Hinheiten sind: h =c =1; Nukleonenmasse M = 1. Dann gibt es nur 
die Dimensionen [Lange] und [Linge]-1. 

Die Fermische Idee war, da8 beim ZusammenstoB etwa zweier hochenergetischer | 
Protonen sich in sehr kurzer Zeit (etwa 1/u, wenn u die Masse des z-Mesons ist) 
und in sehr kleinem Volumen (etwa (42/3) (1/)? ~Nukleonvolumen) ein Gleich- 
gewicht zwischen den méglichen Endzustanden einstellt. Wenngleich im ersten 
Augenblick der Kollision nur Materie (ohne weitere Struktur) in dem kleinen 
Volumen existiert, so soll doch nach der Zeit 1/u eine Struktur sich gebildet 
haben, die es erlaubt, von Teilechenzahlen und Impulsen zu sprechen. Wenn 
das kleine Volumen, in dem die Wechselwirkung stattgefunden hat, verlassen 
BS ist diese Struktur vollstindig festgelegt. Betrachten wir eine bestimmte 

eaktion: 


p+p-N+N-4 nz. (1) 
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Das heibt: Beim ZusammenstoB zweier Protonen entstehen zwei Nukleonen und 
nm Pionen. Die Ladungsverteilung soll uns im Augenblick nichts angehen. Im 
Endzustand ist n festgelegt, ebenso sind es die Impulse aller Teilchen. Aber bei 
mehrfacher Durchfiihrung desselben Experiments werden diese Daten immer 
wieder anders ausfallen. Fermr geht nun davon aus, daB (ahnlich wie in der 
statistischen Mechanik) die Wahrscheinlichkeit eines Endzustandes proportional 
ist der Anzahl der Méglichkeiten, ihn zu verwirklichen. Diese ist gegeben durch 
die Anzahl der quantenmechanischen Zustinde, die dem gleichen Endzustand 
angehoren. 

(Im Sinne der Quantenmechanik ist das, was hier Endzustand genannt wird, 
noch eine Zusammenfassung von vielen quantenmechanischen Zustianden). 
SchlieBlich ist, wenn man die Teilchen durch ebene Wellen beschreibt, die Anzahl 
der Zustainde in einem Normierungsvolumen V gegeben durch die Anzahl von 
ebenen Wellen, die die Randbedingungen erfiillen. Bekanntlich erhalt man diese 
Zahl auf einfache Weise, indem man das ganze Phasenraumvolumen (zwischen 
den Energieschalen Hy und L,+ AE; E, ist die festgehaltene Gesamtenergie) 


yare. fil. f dp, ...b Pats 
fig HAE 


durch die GréBe der Zellen h3"+?) = (27)3™+2) (da & = 1) dividiert. 

Da der Endzustand bereits festliegt, sobald die Teilchen ein Volumen der GréBen- 
ordnung (42/3) (1/u)? verlassen, sollte man V ~ (42/3) (1/u)? setzent). Fermi 
argumentiert weiter, daf die Impulserhaltung die Zahl der Freiheitsgrade um 
eins vermindert und gelangt so zur Aussage 


1 3(n+ 1) 4zn 1 m+ 
P,~ (5. (5 a) | B= Sa) dp. dPaey (2) 


250 


wo P,, die Wahrscheinlichkeit fiir das Entstehen von n Pionen und ¢; = V p} + mi 
die Energie des i-ten Teilchens ist. Im Abschnitt B dieser Arbeit soll gezeigt 
werden, daf man zu einer prinzipiell ahnlichen Formel mit Hilfe der S-Matrix- 
theorie kommen kann. Dabei ergeben sich von selbst Verfeinerungen und eine 
vertiefte Einsicht. Wir wollen daher hier nicht versuchen, die intuitive Begriin- 
dung, die soeben gegeben wurde, zu verbessern. Stattdessen sollen einige Ur- 
sachen (und ihre Beseitigung) dafiir erwahnt werden, daf} (2) sich nicht sonderlich 
gut bewahrt hat. Diese Ursachen sind von mehreren Autoren beseitigt worden 
und es hat immerhin einige Jahre gedauert, bevor die verbesserte Theorie erfolg- 
reich war: 

4. Es ist méglich, da auch nach Verlassen des Wechselwirkungsvolumens 
(42/3) (1/u)? noch gewisse Wechselwirkungen weiterbestehen: Resonanz- 
zustande oder mit anderen Worten auferst kurzlebige instabile ,,Teilchen™, 
die durch starke Wechselwirkung in Zeiten zerfallen, welche nur wenig langer 
sind als die ,,Elementarzeit“ 1/u. Diese ,,Teilchen® spielen aber bei der Ein- 
stellung des ,,Gleichgewichts“ voriibergehend die Rolle eines echten Teilchens 
und sind als ein solches zu behandeln [3, 4, 5, 6]. 


1) Dieses Argument stammt nicht von Fermi und wird hier nur zur Plausibilitatsbetrachtung 
benutzt. Spater wird es durch ein besseres ersetzt werden. 
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. Die oben erwihnten sehr instabilen Teilchen sind nicht beobachtbar. Man 
muB ihren Zerfall beriicksichtigen und die Zerfallsprodukte (und deren Spek- 
tren) zum Endzustand hinzufiigen. Das gleiche gilt fiir andere instabile Teil- 
chen, die auf Grund elektromagnetischer (x°—> 2 etc) oder schwacher Wechsel- 
wirkungen zerfallen, falls sie — wie das x® — iiberhaupt nicht beobachtbar 
sind oder falls man sich gerade fiir die Zerfallsprodukte (etwa Neutrinos) 
interessiert [3, 6, 7]. 

. Wenn Teilchen-Antiteilchenpaare gebildet werden, kénnen solche Paare das 
Wechselwirkungszentrum mit so kleinem Relativimpuls verlassen, daf sie 
lange beieinander bleiben und sich spater gegenseitig vernichten. Die Wahr- 
scheinlichkeit dafiir ist berechenbar [8] und der beobachtbare Endzustand 
kann entsprechend korrigiert werden. 

. Die Erhaltung des isobaren Spins und die Ununterscheidbarkeit gleicher 
Teilchen mu streng beriicksichtigt werden. Dabei treten Faktoren auf, 
die leicht sehr groB werden und die relativen Verhaltnisse der Wahrschein- 
lichkeiten vollstandig verandern kénnen. Unter Annahme der Ladungsunab- 
hangigkeit kann man diese Faktoren berechnen [9, 10, 19]. 

. Das in (2) auftretende Integral bekommt, wenn man die Impulserhaltung nicht 
durch Weglassen einer Integration, sondern streng beriicksichtigt, fiir n Teil- 
chen die allgemeine Form 


eX (E, my, ... mq, P) = { 6(E — S &:) 6(P — Spi) dp, ...d py. 


Hs 1a8t sich nur fiir n < 3 streng ausfiihren [//]. Fiir » > 3 kann man es 
fiir jeden der beiden Grenzfiille 


a) Mm > Dis ei = Vp} +m ~m,+ = : (unrelativistisch) 


am; 
b) m<M3.e= V pi +m ~ p; (ultrarelativistisch) 


exakt lésen. Hat man ein System mit schweren (Nukleonen) und leichten 
(Mesonen) Teilchen, so liegt die Versuchung nahe, fiir die schweren Teilchen 
Fall a) und fiir die leichten b) anzunehmen und ferner zu glauben, daB jedes 
System fiir sich (im Schwerpunktsystem aller Teilchen) Gesamtimpuls Null 
habe und die Totalenergie sich so auf beide Teilchengruppen verteile, daB 
L,+ B,= EF und Fyn = Lying ist. Man hatte dann etwa (IM heiBt: groBe 
Massen) 
Q* ou (£,,P =0) x Om=0 (£,,P = 0), 


und in der Tat ist diese Naherung, die grobfalsch ist, gelegentlich benutzt 
worden. Besser wire es schon (bei der gleichen Naherung fiire = V2 a m?) 
wenigstens nicht anzunehmen, daf die beiden Teilsysteme Energie und Impuls 
in so spezieller Weise untereinander aufteilen. Dann hatte man iiber alle mog- 
lichen Verteilungen zu integrieren und erhielte 


ot = [dH dP’ oh (B — B', P’—P) - oho (EP). 


Es konnte gezeigt werden [12], da auch diese vielbenutzte Naherung leicht 
vollkommen falsch werden kann; die Abweichungen vom exakten Wert 
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kénnen ohne Weiteres einen Faktor Fiinf erreichen. Das liegt daran, daB sich 
in Wirklichkeit die Teilchenimpulse im Mittel so einstellen, daB weder a) noch 
b) eine gute Naherung ergeben. Ferner wirkt sich ungiinstig aus, da® nur die 
schweren Teilchen mit ihren Massen explizite eingehen, wihrend die leichten 
alle die Masse Null erhalten. Zu welcher Gruppe man die K-Mesonen zihlen 
soll, bleibt fraglich. 

Wenn man nun schon intermediire instabile Teilchen einfiihren und auch die 
Produktion von K-Mesonen und Hyperonen berechnen will, so muB man, da 
sich beziiglich des Phasenraumintegrals 9* diese Teilchen nur durch ihre 
Massen unterscheiden, schon zu besseren Methoden greifen und o* exakt oder 
wenigstens so genau berechnen, daf} die Massenunterschiede richtig beriick- 
sichtigt werden. Dies wurde durch eine Sattelpunktmethode [13] und durch 
eine Reihenentwicklung [3] versucht. Eine Methode, die weiter unten aus- 
fiihrlicher beschrieben wird, benutzt eine elektronische Rechenmaschine und 


ein Monte-Carlo- Verfahren [12, 14], wobei die exakten Ausdriicke ¢ = Ver+ m? 
benutzt werden. Dies hat gegeniiber den anderen Techniken den Vorteil, da8B 
man bis zu 16 Teilchen im Endzustand haben kann, deren Massen auf 101) 
verschiedene, beliebige Werte verteilt sein kénnen. Gleichzeitig mit o* erhalt 
man die Energiespektren aller Massen und eine Fehlerabschatzung. Man 
laBt die Maschine solange rechnen, bis ein vorgeschriebener Fehler unter- 
schritten ist. Einige Prozent wahrscheinlicher Fehler gegeniiber dem exakten 
Wert werden in wenigen Minuten Rechenzeit erreicht (CERNs Ferranti- 
Mercury Computer). 


6. Da nicht alle ,,starken‘‘ Wechselwirkungen gleich stark sind, ist es zweck- 
maBig, verschiedene Wechselwirkungsvolumen Q,, 2, ... einzufiihren [14, 8]. 


7. Wenn man sich bei gegebener Gesamtenergie alle energetisch méglichen und 
mit Auswahlregeln (Isospinerhaltung, Baryonerhaltung, Strangeness-Er- 
haltung) vertragliche Endzustiinde (definiert durch Teilchenzahlen) auf- 
schreibt, also etwa 

N+N+272 
1 Sil reg Mie a oe 
“~NI+N+N+A+4+K+4+32 
etc., 


so sieht man leicht, daB& die Zahl der méglichen Endzustande mit wachsender 
Gesamtenergie stark zunimmt. Zaéhlt man nur diejenigen, die nennenswert bei- 
tragen, so findet man etwa fiir p-p-StoBe 


10 Endzustande bei 2.75 GeV Primarenergie (Laborsystem) 
50 ey ” 6.2 GeV 29 ”» 
600 Fe 20. GeV 95 rp 


Jeder solche Endzustand enthalt im allgemeinen 3 bis 6 Teilchen verschiedener 
Masse und zu jeder Masse gehért als Ergebnis der Monte-Carlo-Rechnung ein 
Energiespektrum, welches (als Histogramm) etwa 30 Zahlen enthalt. In 
25 GeV Fall hat man rund 600-5 = 3000 Spektren zu je 30 Zahlen. Diese 


1) Die Zahlen 16 und 10 sind nicht durch die Monte-Carlo-Methode bestimmt, sondern durch 
das gewahlte Rechenprogramm. Sie kénnen im Prinzip erhoht werden, aber das ist praktisch 


unndotig. 
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Spektren miissen mit den Wahrscheinlichkeiten der Endzustande, zu denen 
sie gehéren, multipliziert und dann, nach Teilchen geordnet, tiberlagert 
werden, um die Gesamtspektren der verschiedenen Teilchensorten zu ergeben. 
Das ist eine Aufgabe, die nicht am Schreibtisch erledigt werden kann. Haufig 
hat man daher einfach einen Teil der Endzustande fortgelassen und mittlere 
Teilchenzahlen und: Spektren aus wenigen, als typisch angesehenen Endzu- 
stinden berechnet. Dabei kénnen, abgesehen von unzuldssigen Naherungen 
fiir o* und Nichtbeachtung anderer statistischer Faktoren (Isospin, Ununter- 
scheidbarkeit, Spin), weitere betrachtliche Fehler entstehen. Beispiel: Man 

findet fiir 25 GeV p-p-StoBe fiir die mittlere Pionenerzeugung*) 


5.385 wenn nur Endzustinde 2N + nz betrachtet werden, 
4.23 wenn Nukleon-Paarerzeugung hinzugenommen wird, 
3,37 wenn alle beitragenden Endzustande zahlen. 


Entsprechend andern sich die Pion-Spektren. 

Wenn man einigermafen zuverlissige Aussagen haben will — oder besser: 
Wenn man die wirklichen Aussagen der statistischen Theorie haben will, 
mu man alle wesentlich beitragenden Endzustande einschlieBen. Dabei stellt 
sich im allgemeinen erst beim Rechnen heraus, welche Zustande mitzunehmen 
sind. Die oben angegebenen Zahlen sind beim Durchrechnen der betreffenden 
Falle praktisch ermittelt worden. 

Damit steht aber fest, da fiir mehr als 3 GeV Primarenergie die Auswertung 
(d. h. die gewogene Uberlagerung der Spektren) mittels elektronischer Rechen- 
maschinen erfolgen mu’. Ein Programm, welches das automatisch und sehr 
flexibel durchfiihrt [16, 17], wird weiter unten in groBen Ziigen beschrieben. 


8. Es ist wiinschenswert, die statistische Theorie soweit irgend méglich von 

modellmaBigen Vorstellungen zu befreien. Man wird also versuchen, von der 
S-Matrix her eine neue Begriindung zu finden [78]. Dabei miissen sich die 
oben erwahnten weiteren statistischen Faktoren von selbst einstellen und 
das Wechselwirkungsvolumen eine bessere Interpretation erfahren. 
Freilich zeigt sich, da8 Modellvorbestellungen nicht ganz zu vermeiden sind — 
eine triviale Folge der Tatsache, dafs im Gebiet der starken Wechselwirkungen 
keine Feldtheorie bisher in der Lage war, S-Matrixelemente explizite zu be- 
rechnen. 


Im folgenden Abschnitt B soll die Neubegriindung der statistischen Theorie aus- 
gehend von der S-Matrix kurz geschildert werden. Dabei werden die Punkte 
1, 2, 3, 4 und 6 der obigen Aufzihlung automatisch mit erledigt. 

Im Abschnitt C wird einiges iiber die Rechenmethoden gesagt, womit die Punkte 5 
und 7 der obigen Liste ebenfalls behandelt sind. 

Endlich werden einige Anwendungen und im Flu8 befindliche Weiterentwick- 
lungen gestreift (Abschnitt D). 


B. Ableitung der statistischen Theorie aus der S-matrix 
Die folgenden Ausfiihrungen sollen nur einen Uberblick geben. Fiir Einzelheiten 
siehe [18]. 


') Diese Zahlen sind nicht fiir Paarvernichtung im Endzustand korrigiert, sie dienen nur zur 
Illustration. 
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Die Wahrscheinlichkeit, von einem Anfangszustand |i> in irgendeinen einer 
Menge F von Endzustéinden zu gelangen, ist gegeben durch 


PRS SISO Ps (4) 
fer 


wobei ein konstanter Faktor weggelassen wurde. Wenn wir die Endzustinde 
< f| durch die Massen und Impulse sowie Spins und Isospins der herauskommen- 
den Teilchen charakterisieren, so folgt (7';, t; sind Isospin und dessen dritte Kom- 
ponente) 


P(f)= eh Py th---Tatn, Pr---Pny M-..Mn|S|t> PP dp,...dpn. (5) 
; F 


Man kann das quadrierte Matrixelement ,,reduzieren‘‘, indem man von der 
Jsospinerhaltung Gebrauch macht: 


Ber alas Deis Dns M20 %y, |S 2. >|? = 
Soe Faas. Ug, tye tah eo Dy se Das Why oe sty |S’ le a: ceo, 18) 


wo O7(T,...7',) ein wohldefinierter Koeffizient ist. Indem man auf die Kenntnis 
der Ladungsverteilung verzichtet (was aber nicht sein muB, siehe [7, 19]) kann 
man tiber alle méglichen t¢;-Konfigurationen summieren (das ist mit dem Summen- 
zeichen in (5) gemeint) und erhalt die Gesamtzahl von Moglichkeiten, aus n Teil- 
chen mit Isospins 7',...7',, eine Isospineigenfunktion |7’,¢> zu konstruieren, wo 
T', t den Isospin des Anfangszustandes 7 charakterisieren. Da uns hier nur Teilchen 
mit Isospins 7 = 0, 1/2, 1, 3/2 interessieren und die (7 = 0)-Teilchen keinen 
EinfluB haben, nennen wir diese Anzahl von Moglichkeiten [9, 10, 19] 


Wapy(T'); a = Zahl der Teilchen mit Isospin 1/2 
6 = Zahl der Teilchen mit Isospin 1 
y = Zahl der Teilchen mit Isospin 3/2 
T = Isospin des Anfangszustandes. 


Ferner ziehen wir eine 6-Funktion fiir Energie- und Impulserhaltung aus dem 
quadrierten Matrixelement heraus und finden 


P(F) = Wapy (2) {|<Pr---Pn a---1in| 82>? 0 (BS 6) (YP) dPy---dPns (7) 


-wobei wir die relativistische Invarianz der S-Matrix benutzt haben, um ein fiir 
allemal in das Schwerpunktsystem »’ p; = 0 zu gehen. 

Da P(F) und W,,, (2) dimensionslos sind, aber die 6-Funktionen und die 
Differentiale dp,...dp, nicht, muB das verbleibende quadrierte Matrixelement 
die Dimension 4 #-V”~1 haben, wo 4Z# eine Energie und V, ein N ormierungs- 
volumen ist, das spater wieder herausfallen mui, da seine Grobe irrelevant ist. 
Diese Faktoren ziehen wir heraus, um wieder ein dimensionsloses quadriertes 


Matrixelement zu haben: 


P(P) = Wapy(T)-AE- V2-*{ | <py..-| 8" Pd (2 — S'ei) (SP) EPr--- Pn. (8) 
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Jetzt beriicksichtigen wir noch die méglichen Spineinstellungen durch den Faktor | 


nN 


IT (24; +1), 


j=l 


wo J; der Spin des j-ten Teilchens ist1). Ferner miissen wir die Ununterscheid- 
barkeit gleicher Teilchen in Rechnung stellen. Beim Integrieren kommen alle 
Permutationen vor und infolgedessen werden alle Zustéinde mehrfach gezahlt. 
Das kompersieren wir durch einen Faktor 


1 
TN!’ 
a 


wo N,...N;, die Anzahlen der Teilchen der Sorten 1... sind ohne Ansehen 
von Spin und Ladung (etwa NV, = Zahl aller Pionen). 


a 


A V nd 
Ee em Hee Wear hye lane -AEH x 


x [py ...[8" |i)? 6(B— Zei) 6(Zp,) dp, ... apn (9) 


Das ergibt {wir setzen Vy = 


Wir setzen zur Abkiirzung voriibergehend fiir den Faktor vor dem Integral 


IT (24; + 1) V n—1 
= + Wasy ( )- | aes AE 


und nehmen an, wir wollten das Spektrum des n-ten Teilchens berechnen. Man 
erhalt es offenbar durch Nichtausfiihrung der betreffenden Integration, (dp,= 
= 27p,2 dp, sin 3, + d%,, die Polarachse sei durch das ankommende Primar- 
teilchen gegeben): 


P(F, pp) dp, = 2xC'p; sin 9, dp, 9, { |(Py --. Pal S’” |4)[2.5(E — Lei) x 
x 6 (2'p;) dp, sei Apn—1 O (11) 


Da wir iiber das quadrierte Matrixelement so gut wie nichts wissen, ziehen wir 
seinen Mittelwert vor das Integral: 


P(F, pp) dp, = 2nC py sin B, dp, > dO, > S(E, T, ny... Ty... Tas Pas On) X 
x [ 6(H — e, — X’e) 6(p, + Z'p,) dp, ... dpy—y (12) 


Die Funktion S ist der erwaihnte Mittelwert, 7,, 7 ... sind alle anderen Parameter, 
die den Mittelungsprozess iiberlebt haben. Wenn wir § als Funktion aller Varia- 


1) Warum Spin und Isospin verschieden behandelt werden, ist in [18] ausfiihrlich diskutiert. 
Der im allgemeinen groSe Bahndrehimpuls ist die Ursache. 


(10) 


| 
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blen berechnen kénnten, wiren wir in der Lage, die Wahrscheinlichkeit dafiir an- 
zugeben, da genau die Teilchen 1, 2 ... » (mit Massen My --. My, Spin und Isospin 
J,...dy baw. T, ...7',) entstehen und daB das Teilchen mit der Nummer n einen 
Impuls zwischen p, und p, + dp, hat und einen Winkel zwischen O,und 3,-+ dd, 
mit der Richtung des Primirteilchens bildet. 

Es kann keine Rede davon sein, da man die Funktion S so weitgehend berechnen 
k6énnte. In der Tat ist (12) eine Aufspaltung des Problems in Faktoren, von denen 
die Funktion S das Wesentliche iiber die Wechselwirkung enthiilt, wihrend alle 
ubrigen eigentlich trivial sind. Da8 die trivialen Faktoren aber den numerischen 
Wert von P(F) weitgehend bestimmen und iiberdies nicht leicht zu berechnen 
sind, werden wir noch sehen. 

Da man noch keine Méglichkeit hat, selbst fiir kleine Teilchenzahlen — geschweige 
fiir groBe — S-matrixelemente explizit auszurechnen, mu8 man im folgenden 
bescheidener sein. 

Zunachst wollen wir auf die Frage nach der Winkelverteilung verzichten, indem 
wir iiber den Winkel #, integrieren. Das hat zur Folge, daB kein Widerspruch 
zwischen den Aussagen der Theorie und einer beobachteten Winkelverteilung 
im Schwerpunktsystem bestehen kann. Ein oft gehértes Argument gegen die 
statistische Theorie wird damit entkraftet: Eine isotrope Winkelverteilung wird 
in der hier gegebenen Begriindung der statistischen Theorie weder vorausgesetzt 
noch gefolgert. 

Zur Durchfiihrung der Winkelintegration definieren wir 


o*(H, m, ... My, P) = [ 6(H — Ze) (Zp;— P) dp, ... dp, (13) 


=i 
leah -o* ist in unseren Einheiten die Zahl der Zustaénde im Kasten vt 
MLA 


Da iiber alle p; integriert wird, kann o0* offenbar von der Richtung von P nicht 
abhangen (also auch nicht vom Vorzeichen). Schreiben wir also (12) in der Form 


P(F, pn) dp, = 2x Op2 dp, sin 0, dd, + S(E, T|my ... tMtq, «+ Dns On) X 
x o* (L — En, My +--+ Mn-13 Pn)» (12) 


so bleibt o* von der Winkelintegration unberiihrt und wir erhalten [mit etwas 
yveranderter Bedeutung von P(F, p,)] 


Php.) ap, = CSL, Titty Noe noes Us 00-2 ins Da) 
x oF (E — Eq, .-- sPn) 40 Dn ADy- (12’’) 


Wir behaupten nun, da® dieser Ausdruck (der noch exakt richtig ist) mit Hilfe 
plausibler Naherungen fiir S ziemlich gut berechnet werden kann. Dabei nehmen 
wir fiir den Augenblick an, C und o0* seien gegeben. Wir kiimmern uns nur um 
die Funktion S und diskutieren die méglichen Naherungen: 
«) Der Faktor ph: o*(E — &, m1... ™Mn-1, Pn) wird zu Null fir p, = 0 und 
0 GRE Se is i in Maximum. 
fir ¢, =Von tm =F — (m+ .---Ma-1); dazwischen liegt ein 
Dieser walle tragt der Kinematik (Energie-Impulsbilanz) und der Kin- 
schrankung der Freiheitsgrade der (n — 1) restlichen Teilchen durch das n-te 
Rechnung. S enthalt also in seiner p,-Abhangigkeit diese Effekte nicht mehr. 
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Daher scheint es plausibel, daB S nur noch schwach von p, abhangt und daher 
im Vergleich zum stark verdnderlichen p,-e* als Konstante anzusehen ist. 


Wir streichen also p, aus S. 

6) Der Faktor C [siehe (10)] zahlt unter anderem die Méglichkeiten, mit n 
Teilchen der Isospins 7... 7, und beliebigen t,...t, emen Gesamtisospin 
T, t za kombinieren. Die Zahl dieser Moglichkeiten variiert stark und wir 
werden annehmen, da demgegeniiber S von 7, ... 7, nur so schwach abhangt, 
da& diese Variablen ebenfalls zu streichen sind. Die Abhangigkeit von 7’ zu 
streichen scheint weniger geraten. Entweder man beschrankt sich aut An-- 
fangszustande, die Isospineigenfunktionen sind (etwa pp) — dann gehoren alle 
Endzustaénde zum gleichen 7 und ihre relativen Gewichte werden weitgehend 
T-unabhangig. Oder man hat gemischte Anfangszustinde (pn, 7 p,pp) — dann 
muB8 man Zusatzannahmen machen. 

Wir beschrénken uns auf reine Isospin-Anfangszusténde und streichen 7’ aus 


S. 
y) Jetzt lautet (12’’): 


PF) dpa = CS (Lyng he -.-) 40 Pah Dp 408 i bn Ma ne 


Um eine Abschitzung fiir S zu bekommen, integrieren wir zunachst iiber p, (was 
an § nichts andert). Aus (13) folgt 


[Aap o* (EE — &4, My... Mq_1) Pn)EDn = 0* (E, My... Mp; 0) (13’) J 


[ P(P)dp, =P (iy ng») = C8o* (Bm No ---, 0) (14) 


ist jetzt die Wahrscheinlichkeit, genau n Teilchen mit den gegebenen Eigenschaften 
(Massen, Isospin, ete.) aus dem festgelegten Anfangszustand (H, 7’) zu erzeugen. 
Diese Wahrscheinlichkeit ist gewonnen aus der Vorschrift: Man integriere das 
quadrierte S-matrixelement tiber alle zum gewiinschten Endzustand beitragenden 
differentiellen quantenmechanischen Zustinde (Impulse, Ladungen, etc.) und 
diese Wahrscheinlichkeit ist in (14) geschrieben als: [Mittelwert des quadrierten 
Matrixelemente] mal [gesamtes Integrations (Summations) Intervall], wobei S 
der Mittelwert ist. 

Nun betrachten wir den inversen Proze8: die n Teilchen des bisherigen End- 
zustandes sind in einem Kasten (Normierungsvolumen V) enthalten und wir 
fragen nach der Wahrscheinlichkeit dafiir, daB® sie sich alle wahrend einer Zeit 
der GroBenordnung 1/u in einem Volumen der GréBenordnung Q, = (42/3) (1/u3) 
treffen und miteinander zu zwei Protonen verschmelzen, die mit hoher Energie 
das Gebiet der Wechselwirkung verlassen. 

Da fiir starke Wechselwirkungen die Invarianz unter Zeitumkehr angenommen 
werden darf, ist das quadrierte S-matrixelement dasselbe wie beim zuerst be- 
trachteten Prozef. Der jetzige Endzustand (d.i. der frithere Anfangszustand) 
ist ein reiner Zustand, daher entfallt die Integration iiber Endzustande. Der 
jetzige Anfangszustand (d. i. der friihere Endzustand) ist aber ein Gemisch aus 
vielen differentiellen Zusténden (Impulse, Ladungen, etc.) tiber die man das 
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quadrierte Matrixelement zu mitteln hat. Nach unseren obigen Uberlegungen 


ist aber die Funktion S(H, 7,7...) genau dieser Mittelwert. Also ist 


Pinvers = S(E, Unt Noes). (15) 


Elementare Statistik lehrt, daB diese Wahrscheinlichkeit proportional sein muB 
zu ({2)/V)"~1. Wirschreiben n — 1, weil der Ort des Zusammentreffens innerhalb V 
gleichgiiltig ist, also nur n — 1 Teilchen ein beliebiges treffen miissen. 

Fir die Energieabhingigkeit kann man folgende Annahme machen: Damit 
zwei Protonen von hoher Energie herauskommen, die im Schwerpunktsystem 
durch Lorentzkontraktion abgeplattet sind (ihr Radius in Flugrichtung ist 
I/yu, wo y =2/E und E die totale Schwerpunktenergie ist), mu8 im letzten 
Augenblick vor ihrer ,,Geburt“ auch das Wechselwirkungsvolumen schon ab- 
geplattet sein, da es aus zwei sich mit hoher Geschwindigkeit lésenden Teilvolumen 
besteht. Also miissen die n Teilchen sich im Volumen Q treffen: 


CAF 


2 3 
are 2Q,=—>— (—-] ; # = totale Schwerpunktsenergie. 


ic 


6) Endlich kann man die Abhangigkeit von den Parametern 7, 7... wenigstens 
andeutungsweise beriicksichtigen. Die Massen z. B. sind durch Wechselwirkung 
hervorgerufen, und verschiedenen Massen entspricht die Existenz verschiedener 
Wechselwirkungen mit verschiedenen Kopplungskonstanten und Reichweiten. Es 
liegt daher nahe, den Ausdruck (2/Q)” zu ersetzen 


aN (= ; F 
(2) -(2).(29.. 

wo zwar die GréSenordnung aller 2; ungefahr festliegt, aber doch Unterschiede 
bis etwa zu Faktoren 10 entsprechend der jeweiligen Wechselwirkung vorkommen 
k6énnten. Die obige Begriindung legt es nahe, 2 als Ausdruck der Reichweite 
der Wechselwirkung anzusehen. Die Erfahrung zeigt aber, dafi wenn man die 
Q; fiir alle vorkommenden Massen durch (42/3) (1/m,)? - 2/H definiert, keine Uber- 
einstimmung mit dem Experiment erzielt werden kann. Wenn man an die St6- 
rungstheorie denkt, so wiirde man fiir jedes erzeugte Teilchen einen Faktor gj 
(Kopplungskonstante) erwarten und dann hiatten Nukleonen, Antinukleonen 
und Pionen das gleiche 2. Es erscheint daher zweckmabig, die 2; als Parameter 
anzusehen. Wir werden, um nicht zu viele Parameter zu haben, drei verschiedene 
Q zualassen und dem Q,, den in (16) gewahlten Wert geben: 


8 oO n—1 LOP Ny (Oe Ny : Qs; oe 
V Aa V Palme 


dann kénnen wir die beiden anderen fiir Hyperonen, K-mesonen, Antinukleonen, 
Antihyperonen, gebundene Zusténde wie Deuteronen, H® He® oder Antideutero- 
nen usw. und eventuell Resonanzzustiinde (Isobaren) verwenden. 

Da Pinvers ZU (2/V)"—1 proportional ist, setzen wir 


443 O Ny O Ne o Ns 
S(£, 4 Yor 5) = A v (7) Y (7) Z (7) 5) (16) 
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wo A noch eine — wie wir hoffen — langsam veranderliche Funktion der Energie 
E, der Massen und aller anderen unterdriickten Variablen ist. Wir werden A als 
Konstante behandeln. 

Tragen wir (16) und (10) in (12’’) ein, so folgt 


LTT ae Ny Ds Ne Q; Ns 
Pur pbtes = pre Woon aes] “Leal [eas] * 
k 


x o* (EH — &j, My... Hj. My, Di) 4p; Ap; (17) 


(22) 22)8 22)8 


wo wir statt des n-ten Teilchens jetzt allgemeiner ein beliebiges, 7, herausgegriffen 
haben. #; heiBt: die Masse m, fehlt. ' 
In (17) sind alle GréBen wohldefiniert und berechenbar. Die Konstanten A und 
AE lassen wir fort, da die Wahrscheinlichkeiten ohnehin durch Division mit 
> P(F) zu normieren sind. 

F 


Wir bemerken nun, daB instabile Teilchen, deren Lebensdauer zu kurz ist, als 
da8 man sie beobachten kénnte, aber doch hinreichend lang, um das Gebiet 
der Wechselwirkung vor ihrem Zerfall zu verlassen, in (17) ganz zwanglos als 
Teilchen eingefiihrt werden kénnen. Nehmen wir als Beispiel das Nukleon-isobar - 
(x-N Resonanz 7’ =J = 3/2; m*=~ 1.31). Man fihrt es als ein Teilchen der 
Masse m* ein und berechnet sein Spektrum mit Hilfe von (17). 

Man konnte fiir jedes derartige ,,Teilchen“ ein spezielles Wechselwirkungsvolumen 
Q* ansetzen, denn der inverse ProzeB kann so vorgestellt werden: 


N YE 


Fiir die Vereinigung x -+ N-—-+N* schreibt man Q* und fiir die Vereinigung 
N+ N*+27-+...+ 7 ein weiteres Produkt von Q’s. Wenn man nun QF Q; 
wahlt, kommt man zu dem von LINDENBAUM und STERNHEIMER vorgeschlagenen 
Isobarenmodell [6], in dem Pionen nur auf dem Weg iiber das Nukleon-isobar 
erzeugt werden. Wir wollen das nicht tun, denn bei unserer Deutung von Q 
liegt dazu kein Anla& vor. DaB eine Tendenz zur Bildung eines ,,Isobars‘“‘ (auch 
héherer Masse) besteht, wird schon von den iibrigen Faktoren bewirkt. 

(17) beschreibt das Energiespektrum des i-ten Teilchens im Schwerpunktsystem. 
Geben wir dem Nukleon-isobar (oder irgendeinem anderen instabilen Teilchen) 
die Nummer ?#, so haben wir in (17) das Spektrum der ,,Isobaren‘‘. Da die Lebens- 


dauer von der GréBenordnung 10-22 sec ist, konnen nur seine Zerfallsprodukte be- 
obachtet werden. 
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Nimmt man an, die ,,Isobaren‘‘ mit Masse m* seien im Schwerpunktsystem 
isotrop verteilt (hier muf eine Annahme iiber Winkelverteilung gemacht werden) 
und ihre Zerfallsprodukte mit den Massen m und  seien im Ruhsystem des 
Isobars wiederum isotrop verteilt, so findet man (durch Transformation vom 
Ruhsystem des ,,[sobars“ in das Schwerpunktsystem) fiir das Energiespektrum 
der Massen m bzw. « folgende Form [14] 


VW 
m* * 
Gra oe 
Le m 
s y y 


m* 


E Ey? 
Tae ie fe (= ) S3(mt a i] —m*; B=e+p. (18) 
Hier bedeuten 


w*(y) das Spektrum der kinetischen Energie y des instabilen Teilchens der 
Masse m* (welches im ProzeB m*->m-+ mw zerfallt) im Schwerpunkt- 
system; 

wW,(e) das Spektrum der kinetischen Energie des Zerfallsprodukts der Masse yu 
im Schwerpunktsystem ; 


Eu die totale Energie des Teilchens der Masse ~ im Ruhsystem des Teilchens 
der Masse m* (also des Isobars); 
k den Impuls von m (und w~) im Ruhsystem von m*, 


Offenbar mu man w,(e) und w,,(e) zu den anderen Spektren der Teilchen yu 
bzw. m hinzufiigen. (Man erhalt w,,(e), indem man in (18) iiberall m mit yu ver- 
tauscht). Formel (18) ist natiirlich fiir jeden 2-K6rperzerfall brauchbar und dient 
infolgedessen auch zur Berechnung des y-Spektrums in 7°-> 2y und X°-> A®+ y. 
Als weitere Endzustand-Wechselwirkung kommt die Vernichtung eines zunachst 
gebildeten Paares N+ N oder © +2 oder N+ A usw. in Frage. Mit Hilfe 
von o* kann man ausrechnen, wie wahrscheinlich es ist, daB ein solches Paar 
einen Relativimpuls < 7, hat?), der so gewahlt ist, daB beide Partner der méglichen 
Vernichtung etwa fiir eine Zeit ¢>1/u einen Abstand d <1/u behalten, 
so daB die Vernichtung sehr wahrscheinlich wird. Nimmt man an, daf ein ge- 
wisser Bruchteil [etwa o Vernichtung/ctotal (P;1)] solcher Paare verschwindet und 
die entsprechenden Mesonen dafiir auftreten, so kann man die berechneten Multi- 
plizitaten korrigieren. Die Durchfiihrung einer solehen Rechnung ist in [8] im 
einzelnen dargestellt. Man findet z. B., daB bei p-p-StdBen mit 6.2 GeV Primar- 
energie (Laborsystem) etwa 75% und bei 25 GeV etwa 25% der gebildeten Paare 
sich im Endzustand wieder vernichten und daher nicht beobachtet werden. 
Wir kehren zur Formel (17) zuriick und stellen zusammen, was man mit ihrer 
Hilfe berechnen kann. Bezeichnen wir mit 
F eine Reaktion, charakterisiert durch Teilchenzahlen 
(etwa P+ P+>N+2+4+K-+ 5x), mit 
4 die Teilchensorte (z.B. i = 1, 2,3... entsprechen (19) 
, &, 4, 2... 2, y) und mit 
nF) die Zahl der Teilchen der Sorte ¢ (etwa Pionen) in der Reaktion F (fiir das 
obige F ist nz = 5), 
1) Indem man namlich dieses Paar als ein ,,Teilchen‘‘ der Masse (m,-+ mg) -+ pre1/2 (74 + 7p) 
in das Phasenraumintegral aufnimmt. Beweis in [8]. 
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so stellt (17) das (nicht normierte) Impulsspektrum der Teilchensorte a in der 
Reaktion F' dar und die (ebenfalls nicht normierte) Wahrscheinlichkeit fiir die 
Reaktion F folgt durch Integration tiber p,; 


I] (245+ 1) 


P(F) = +__—_ TIN! 
k 


Waperil)- I eae oF (H, m,.-- Mn, p =90). (20) 


Dann ist 


“=w,’ (p) (21) 


das auf eins normierte Impulsspektrum des Teilchens 7 in der Reaktion F. 


SP(F, p;) ne” SPF): ws” (p) + nf” 


F 
LE i ee SP (22) 
. F F 


das totale, auf mittlere Teilchenzahl normierte Impulsspektrum der Teilchen- 
sorte 7 aus allen Reaktionen, 


LP (F) ni” 
ij Wilde = ei =) (23) 


die mittlere Zahl von Teilchen der Sorte 7, gemittelt iiber alle Reaktionen, 
also die mittlere Multiplizitaét. SchlieBlich kann man noch mittlere Energien, 
etwa 
(exinyi = [ Walp) - e(p)dp|(ni) (24) 

berechnen. 
Die >’ sind iiber alle Reaktionen zu erstrecken, die merklich beitragen; das MaB 

F 
dafiir sind die P(F). 
Bezeichnen wir noch mit 


© =(F,, Fo, ... Fy) (25) 


eine beliebig definierte Menge von Reaktionen (bei gleichem Anfangszustand), 
zum Beispiel alle Reaktionen in denen Antihyperonen vorkommen, so folgt 


(26) 


fiir das Verzweigungsverhaltnis zwischen zwei solchen Mengen von Reaktionen. 
Natiirlich kann man auch in (22), (23), (24) die Summation nur iiber F e@ erstrek- 
ken und hat dann die dort definierten GréBen bezogen auf eine ausgewahlte Menge 
von Reaktionen. 


Neue Methoden in der Statistischen Theorie der Teilchenerzeugung 15 


C. Rechenmethoden 


Abgesehen von den Summationen iiber alle Reaktionen F, die schon bei iiber 
3 GeV Primdarenergie (Lab.) auf Rechenautomaten durchgefiihrt werden miissen, 
gibt es nur noch ein Problem bei der Berechnung, namlich die Ausdriicke 


o*(E, Myz++. Mp, P = 0) und 
oF (H — &;, M1... Wty... May Di) 


die in (20) baw. (17) auftreten. 


1.Berechnung der Phasenraumintegrale mit Hilfe einer Monte-Carlo- 
Methode 


Wir geben hier nur das Prinzip, Einzelheiten findet man in [12]. Wie schon er- 
wahnt, kann das Integral 


o* (E, m,... Mp, P) = [oz — > Vr? + m3) 6(> pi— p) dp,...dp, (27) 


nur fiir n < 3 analytisch ausgefiihrt werden. Wollte man nun eine Art Simpson- 
Forme! anwenden, so fande man, da infolge der unendlichen Integrationsgrenzen 
und der hohen Dimensionszahl nur mit groBem Rechenaufwand einige Genauigkeit 
zu erzielen ist. Man kann die Zahl der Integrationen auf zwei reduzieren, tauscht 
dafiir aber sehr unbequeme Integranden ein und muBf singulaére Funktionen 
numerisch integrieren. Unter diesen Umstanden scheint eine Monte-Carlo-Methode 
angebracht, die wir jetzt vorbereiten: 

In (27) trennen wir die Winkelintegrationen ab, indem wir setzen 


Pi=Viei; lei) =1 (28) 
dp, = pidp,de; = D; sin 6; d0; dp; dq, 
dann wird 
Oo (LH, m,...M_,p) = f Pol dp,...dp, p?... Pr(E— > yp + m) x 
0 0 
(29) 
x f...fo(P — XY pei) de... den. 
Richtungen 
Die Funktion 
il n 
w, (P, Dy <-< Pay = (a) -{...f6(P — 3 me) de... de, (30) 
4a Richtungen 


hat eine einfache Bedeutung: w,(P,p,--.Pn)dP ist die Wahrscheinlichkeit 
dafiir, daB n Vektoren gegebener Lingen 7,...p, aber statistisch verteilter 
Richtungen e,...e, sich zu einem resultierenden Vektor ae p;, zusammen- 
setzen, der in dP liegt. Diese Bedeutung abt sich aus (30) unmittelbar ablesen, 
der Faktor (1/4 x)” dient zur Normierung, so daf 


fw (P, Pie+ Pn) EP = 4m [ P2w,(P, Pr... Pr) dP = 1. (31) 
0 
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Aus der anschaulichen Bedeutung von w, ergibt sich ohne Rechnung 
Wy (0, Dye++ Pn) = Waa (Pj Pies» Bie» Pn)s G=1L-...D, (32) 


so daB es geniigt, w,(0,p,... pn) zu berechnen. 
Durch Einfiihren der Fourierdarstellungen der 6-Funktionen in (30) kann man w, 
analytisch ausrechnen. Man findet [12] [sg(x) = Vorzeichen von 2] 


1 


Mia\O> Biss; Pa) 2 (eee rts gS 
1 
a SS On ° Og--- On ° 8 0; ek o;p;)""* |, 33 
Pa Pa Pa [ aan B (2 Pi) (2 o BS 


wobei o,=+1 gt o, = 1) und iiber alle 2-1 méglichen Vorzeichenkombi- 
nationen 05, 03, ... G, Zu summieren ist — eine schéne Ubungsaufgabe fiir das 
Programieren von Rechenautomaten (natiirlich soll w, in méglichst kurzer Zeit 
ausgerechnet werden!). 

Wir sehen w,(P,p,... Pn) als gegeben an. 

Durch Einfiihren der Energien und neuer Abkiirzungen . 


= Vp? + mi > P= Vez — mi 
Ul bs eeeten) = Wall src Da) (34) 
x ( igen 
0 E< Mm 
erhalten wir im Schwerpunktsystem (p = 0) aus (29) 


Q* (2,2... 1My,0) = (4 ae)" [ dey... f dente (2) ... ttn (En) Un (0, £4--. &n) - O(E — SF a). 
My, Mn (35) 


Jetzt gehen wir zu kinetischen Energien iiber: 


und erhalten 


| o* (Hi, m1... Mp, 0) = (4n)" [dty...fdt, ® (t, .. ,) 6(L' — St) mit 
6 0 


| (37) 
D (ty... ty) = Uy (t+ my)... Un(tn + My) - Un (0,6, +'my, ... ty my). 
Zuletzt stellen wir die ¢; als Differenzen von Gré8en 7’; dar: 

tis fe 

(38) 


ee Eh. 
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Damit wird aus (37) 


o* (E,m,... mMn,0) = (aye fav, [aT,... [dt 4 fa, x 
0 TP, Ds Tr—1 


x @(T,,T,— fe eas} 52 6(7— AUR 
Durch Integration iiber 7’, folgt 


D: ae ey 
0% (H, my... my,0) = (4)"[ dT for fare, x 
0 Ts Tn—2 (39) 


SST EU een gee At hs > 27g 


wo ® mit Hilfe von (37) und (34) berechnet werden kann. 

Auf die Formel (39) wenden wir eine Monte-Carlo-Methode an: Wir wahlen 
n — 1im Intervall 0 bis 7’ gleichmaBig (d. h. mit konstanter Wahrscheinlichkeits- 
dichte) verteilte Zufallszahlen und benennen sie der Gréfe nach: : 


1 REE aired REA RT LES 


Einen solchen Satz von n — 1 Zufallszahlen nennen wir eine geordnete Stich- 
probe. Wir behaupten nun: Wenn JN solche geordneten Stichproben 


LA ee ef Ae 4=1...N 


genommen werden, so ist 


ik a i 
faT, [aT,...fdT,4 O(T,,T,— Ty...) = 
0 Ty To 


TO im LY OTO. TH — TO... TO — FO) (40 
— 1 = v O) eo Cee MOD as v 
ae (n nz ieee WH, ( a 08 2 if n ws) ( ) 
Um (40) zu beweisen bemerken wir, dai die Wahrscheinlichkeit dafiir, daB die 
i-te geordnete Stichprobe in die Zellen (siehe Figur) AT,,AT,:.. AT, fallt, 


-—_ a+ a — A ===5 
PAT ALA AL, At vis 


gegeben ist durch 
S95 NI 
(n )! ae : ae = eS = foe ACE AL sae Aly 


'N mal diese Wahrscheinlichkeit gibt die Gesamtzahl der geordneten Stichproben, 
die in genau diese Zellen fallen und alle diese liefern (wenn alle AT';— 0) cen: 


selben Wert von ®. Da ferner 7, < 72... << T ist, wiirden wir die SOE, ®; 
niherungsweise berechnen als = 


N ; hes te ih T 
OT? jew SR far, f ary... [el Ses 
Te O 1 n—-2 


qn—-1 
womit (40) bewiesen ist, da im Limes V — co das Gleichheitszeichen gilt. 
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Damit ist auch prinzipiell klar, wie man auf einem Rechenautomaten g* berechnet: 
Man stellt einen Satz Zufallszahlen her und berechnet ®. Das tut man sehr oft 
und addiert gem&® (40) alles auf. ZweckmaBigerweise laft man N offen und 
berechnet (40) fiir mehrere Satze Sx von je N, Stichproben, die jeweils ein etwas 
anderes 0% ergeben. Der Mittelwert 9* all dieser 9% konvergiert gegen das ,,wahre* 
o* und aus den Abweichungen > (0% — ox)? 14Bt sich der wahrscheinliche Fehler 
K ; 


bestimmen. Diese Probe wird waihrend der Rechnung stets nach jedem Satz von Vg 
Stichproben (mit Ausnahme der ersten paar Satze) durchgefiihrt und die Rechnung 
abgebrochen, sobald ein vorgeschriebener Minimalfehler unterschritten wurde 
[12, 22]. 

Die Spektren aller Teilchen werden hierbei automatisch mitberechnet, indem man 
jedem Teilchen eine Reihe von, sagen wir, 30 Speicherzellen zuteilt, die geeignet 
gewahlten Impulsintervallen A p entsprechen (namlich so, da 30 - 4 p der héchste 
einem einzelnen Teilchen erreichbare Impuls ist). Hat in der k-ten Stichprobe 
das Teilchen i (etwa Pion) den Impulswert p™ = i? + 2tm, erhalten, so 
addiert die Maschine die Funktion (7...) in diejenige der 30 Speicherzellen des 
i-ten Teilchens, in der p liegt (also in die Zelle mit der Nummer [p;/A p]). Be- 
’ trachten wir einmal eine feste Teilchensorte, so sehen wir, daB bei jeder Stich- 
probe die zugehorige Funktion @ in genau eine der fiir das Teilchen reservierten 
30 Zellen Ap addiert wird. Am Ende der Rechnung erhalt man in diesen Zellen 
ein Histogramm dessen Summe iiber alle Zellen (bis auf triviale Faktoren) genau 
o* ist. In der Tat stellt jetzt, wie eine nahere Uberlegung zeigt, dieses Histogramm 
das Spektrum des betreffenden Teilchen dar, d.h. eine Histogramm-Naherung 
fiir 


w;(p) = const- p;? - o* (H — &, my... Wy... Mn; Pi)- (41) 


Diese Naherung enthalt natiirlich in der einzelnen Zelle gréBere statistische 
Fehler als 9*(H,m,...m,,0) selbst (siehe Figur). 

Diese Monte-Carlo-Methode liefert also 

w; (p) in einem Rechengang den Faktor o* (LZ, 

M1... M,,0) und die Spektren w,;(p) 

aller Teilchen. Praktische Rechnungen 

zeigten, daB auf einer schnellen Maschine 

(Ferranti-Mercury bei CERN) Genauig- 

keiten von einigen Prozent in 1—15 

Minuten (je nach Teilchenzahl) erreicht 

werden. Verlangt man, da die Reak- 

tionen mit gro&tem Gewicht (groBes o*) 

iy genau zu berechnen sind (sagen wir 

Hie auf2%), so ergibt sich, da8 die meisten 

anderen (deren 9* viel kleiner ist) nur 

noch auf 10% oder noch viel weniger 

genau bestimmt werden miissen. So konnte man selbst Falle wie 25 GeV (Lab) 

pp-StoBe, zu denen etwa 600 Reaktionen beitragen, in‘ etwa 50—100 Stunden 

Rechenzeit vollstandig erledigen, wihrend bei 6,2GeV nur 8 Stunden _be- 

nétigt wurden, Es wiirde aber aussichtslos sein, etwa > 100 GeV pp-StéBe mit 

diesen Methoden anzugehen. 
Die Brauchbarkeit der Methode scheint wesentlich darauf zu beruhen, daB die 
Zufallszahlen im Mittel zur Gleichverteilung der kinetischen Energie fiihren, 


Ap 
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1.h., fast alle Beitrage zum Integral werden dort berechnet, wo der Integrand 
sowieso ein steiles Maximum hat. Diese Methode entspricht also einer numerischen 
Integration, in der man die Intervallinge dem Verlauf des Integranden anpaBt: 
Kleine Intervalle an steilen Maxima, groBe dort, wo der Integrand klein oder 
langsam. veranderlich ist. 


2. Auswertung und Reduktion der berechneten Spektren und 
Phasenraumintegrale auf physikalische Aussagen 


Die nach der soeben beschriebenen Methode berechneten Spektren und Phasen- 
raumintegrale haben noch keine direkte physikalische Bedeutung. Nehmen wir 
an, die Maschine liefere die auf eins normierten Spektren und die Werte der 
Phasenraumintegrale o* fiir jede Reaktion F: 


of und wi (p) 
sind dann in Form von Lochstreifen gegeben, die fiir die Auswertung wieder in 
die Maschine eingelesen werden kénnen. Sei Lp der Lochstreifen, der of und die 
Spektren w) enthalt. 


Damit die Maschine P(F) fiir alle F berechnen kann, mu8 zusammen mit Lp 
eine Reihe von Daten eingelesen werden, die die nach (20) erforderlichen Infor- 
mationen enthalten, namlich fiir die Reaktion F: 


die Spins J; aller Teilchen, 

die Zahlen Vx gleicher Teilchen, 

a, B, y, also die Zahlen der Teilchen mit Isospins 1/2, 1, 3/2, 
7’, den totalen Isospin, 

Q,, Qs, Q3, die gewahlten Wechselwirkungsvolumen und 

M1, NM, Nz, die zu ihnen gehérenden Teilchenzahlen. 


Ferner miissen fiir jede Teilchensorte (einschlieBlich instabiler Teilchen) wieder 
je 30 Zellen fiir das Spektrum reserviert werden. Das Auswerteprogramm [16, 17] 
sorgt jetzt dafiir, dai das Lochband Ly mit zugehdrigen Interpretationsdaten 
eingelesen, P(F) nach (20) berechnet und alle (auf eins normierten) Spektren 
we ) mit F(P) multipliziert werden, was P(F’, p;) ergibt (siehe[2Z]). Das Band mit 
den Interpretationsdaten enthalt noch die Angaben, an welcher Stelle im Speicher 
P(F, p;) abgelegt werden soll. Das geschieht additiv, nachdem noch mit n(”), der 
Zahl der Teilchen 7 in der Reaktion F', multipliziert wurde. Ferner. wird P(F) 
ausgedruckt und 2P(F) abgespeichert. 

Sind so alle Reaktionen Ly mit zugehérigen Interpretationsdaten eingelesen, 
so befinden sich alle in (22) bis (24) definierten Grofen praktisch fertig zum 
Ausdrucken im Speicher. Man muf nur die Spektren noch normieren und mittlere 
Energien berechnen. SchlieBlich kann man einige der Spektren, die instabilen 
Teilchen entsprechen, ,,zerfallen“‘ lassen, wobei (18) benutzt wird. Wenn das 
geschehen ist, kann man die Spektren noch durch ein geeignetes Interpolations- 
verfahren glitten, um die Reste der statistischen Schwankungen zu beseitigen. 
Da die Spektren (je nach Teilchenmasse) sehr verschiedene Form haben kénnen, 
ist das nicht einfach. Man darf zur Interpolation nicht zu hohe Ordnungen (von 
Polynomen, Fourierreihen oder anderen geeigneten orthogonalen Systemen) 
zulassen, da sonst die statistischen Schwankungen reproduziert werden. Anderer- 
seits kann man Kurvenscharen mit stark wechselnden Formen nicht durch 
Interpolationen niederer Ordnung gut annahern. Daher ist es zweckmafig, durch 


D * 
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geschickt gewahlte Transformationen alle Spektren auf eine Art ,,Hinheitsform“ 
zu bringen, dann zu glaétten und nachher wieder zuriickzutransformieren. Das 
alles mu8 natiirlich automatisch (ohne Hingriff eines menschlichen Beobachters) 
geschehen. Das heift: Das Unterprogramm fiir Interpolation mu die Spektren 
auf ihre Form priifen und individuell jedem Spektrum ,,seine™ Transformation 
zuordnen, die es der gewiinschten ,,Normalform‘ méglichst nahe bringt, mu 
diese Normalform glatten und dann wieder zuriicktransformieren. 

Nachdem all diese Operationen ausgefiihrt sind, 148t man ausdrucken und erhalt 
nur noch eine gut iibersehbare Menge physikalischer Daten. 

Das soeben in grofen Ziigen beschriebene Auswerteprogramm mu8 auf das 
Monte-Carlo-Programm abgestimmt sein. Eine genaue Beschreibung findet sich 
an [ 16; 17 |; 


D. Anwendungen 


Vergleich mit verschiedenen Experimenten 


Die beschriebenen Methoden der statistischen Theorie sind bisher zur Berechnung 
folgender Falle benutzt worden (Energie immer im Laborsystem) 


pp-StoBe bei 2.75 GeV kinetischer Primarenergie [3, 7] 
pp-StoBe bei 6.25 GeV . 5 [8, 29] 
pp-Stobe bei 25 GeV @ ss [18] 
NN-Vernichtung in Ruhe [2/, 22] 

Denteronproduktion in pp-StéBen [31] 


Beim Vergleich der Rechnungen mit Experimenten stellt man folgendes fest: 
1. Bei 2.75 GeV erhalt man fiir das Verhaltnis der Haufigkeiten der Erzeugung 
von 0, 1, 2,3, 4 Pionen gute Ubereinstimmung mit dem Experiment, wenn das 
Nukleon-Isobar (m*= 1.31, 7’ = J = 3/2) als instabiles Teilchen behandelt 
wird; dagegen ist die Ubereinstimmung schlecht, wenn man das Isobar fortlaBt 
[Tabelle 1]: 


Tabelle 1 (entnommen aus [7]): 


Wahrscheinlichkeit in Prozenten 
Endzustand |- 
ohne Isobar mit Isobar Experimentell 
2N 7.13 2.4 unbekannt (klein) 
2N+7 63.9 30.3 36 
2N + 2x 26.9 48.4 48 
2N +32 1.8 18.5 16 
2N+ 42 0.02 0.4 unbekannt (klein) 


Die Tabelle zeigt, wie gro der Hinflu&B des Nukleonisobars ist. Das ist nicht 
selbstverstandlich, denn die Lebensdauer ist von derselben GréBenordnung wie 
die Reaktionszeit: ~10-** sec. Man wird also dieses Isobar bei allen Rechnungen 
beriicksichtigen, sein Wechselwirkungsvolumen sollte dem fiir Nukleonen und 
Pionen etwa gleich sein. 

2a. Bei 6.25 GeV findet man fiir die Produktion von Antiprotonen die richtige 
GroBenordnung (Np): (Nz-) = 10-°, wenn man fiir Antinukleonen dasselbe Wechsel- 
wirkungsvolumen 2 nimmt wie fiir Nukleonen und Pionen. Man kénnte meinen 
dafs an der Antinukleonenproduktion vorzugsweise der ,,core“ beteiligt und daher 


Neue Methoden in der Statistischen Theorie der Teilchenerzeugung 21 


ein entsprechend Kleines 2 zu nehmen sei. Dies fiihrt zum Widerspruch mit dem 
Experiment. Das lé8t vermuten, daf die Nukleon-Antinukleon-Wechselwirkung 
von der Nukleon-Nukleon-Wechselwirkung nicht sehr verschieden ist und der 
beobachtete gréRere totale Wirkungsquerschnitt fiir N-N-Sté8e im wesentlichen 
vom grofien Phasenraum fiir Vernichtung herriihrt. 
2b. Bei 6.2 GeV findet man die richtigen Verhaltnisse fiir die Hiufigkeit von 
Ereignissen mit 2, 4, 6 usw. erzeugten geladenen Teilchen und die richtige totale 
Zahl pro Ereignis [8, 20]. Die Spektren werden aber nicht gut wiedergegeben. 
Durch Einfiihren eines z-z-Resonanzzustandes (,,Pion-Isobar‘ in Analogie zum 
Nukleonisobar), der aus ganz anderen Griinden (Nukleonstruktur, Paarvernich- 
tung) schon gefordert wurde und durch Beriicksichtigung ,nichtzentraler“ 
StéBe konnte man die Ubereinstimmung verbessern [20, 23, 24]. Wir kommen 
darauf zuriick. 
2c. Uber die Produktion von K-Mesonen und -Hyperonen ist nur fiir einige 
Winkel im Laborsystem etwas bekannt. Die Rechnungen in [8] gelten fiir das 
Schwerpunktsystem, eine Transformation auf das Laborsystem wurde nicht 
durchgefiihrt. Es scheint aber, da die Produktion von Strange particles“ 
richtig wiedergegeben wird (Gréfenordnung), wenn man Q, = (1/2 bis 1/5) - Q, 
setzt und fiir die K-Mesonen reserviert, wahrend fiir die Hyperonen 2, zu nehmen 
ist. Mit diesen Werten erhalt man anscheinend auch Ubereinstimmung mit 
Experimenten bei 25 GeV. 
3. Nachdem die Wechselwirkungsvolumen an Experimenten bei niederer Energie 
ungefahr festgelegt sind, kann man bei héheren Energien Voraussagen machen, 
die von Interesse fiir die Planung von Experimenten sind. Aus diesem Grunde 
wurde bei 25 GeV eine Rechnung fiir das CERN Proton-Synchrotron durchgefiihrt. 
Experimentelle Ergebnisse liegen noch kaum vor, da am 25 GeV-Proton-Synchro- 
tron in CERN bisher keine Wasserstofftargets benutzt wurden. Trotzdem hat 
man die Rechnungen (nach Transformation ins Laborsystem) mit Messungen 
verglichen und qualitative Ubereinstimmung gefunden, die teils gut, teils schlecht 
ist. Es scheint aber, da die Abweichungen im wesentlichen damit erklart werden 
k6nnen, daB Eisen oder Aluminium statt Wasserstoff als Targetmaterial benutzt 
wurde [sekundire Reaktionen (Kernkaskaden), Nukleonenbewegung im Kern, 
schlecht definiertes Schwerpunktsystem, usw.]. Die einzige (auf dem Umweg 
iiber Kernemulsionen) an Nukleonen gemessene Zah! ist: Pro Ereignis verlassen 
im Mittel 4.02 (alle St6Be) bzw. 4.68 (pp-StoB) geladene Teilchen die Reaktion 
(private Mitteilung von H. WryzetzEr, Bern), wahrend man fiir pp-StoBe 4.6—4.7 
aus der Theorie [78] erhalt. Vernachlassigt man den Kinflu8 des totalen Isobaren- 
spins, der im pn-Sto8 nicht definiert ist, so lautet die Aussage der Theorie, dai die 
Zahl CO — 2 der neu erzeugten Ladungen etwa 2,6 bis 2,7 ist [siehe Bild 2]. Im pp-StoB 
sollte die Gesamtzahl geladener Teilchen nach dem StoB also etwa 4.6—4.7, fiir 
p-Kern-StoB [= (pn + pp)-StéBe] 4.1—4.2 betragen. Die Unsicherheit des 
theoretischen Wertes kommt von der ungenauen Festlegung von 

Qx = ("/p bis */5) + Qe. 
Interessant wire die Antinukleon- und Hyperonproduktion. Man findet in der 
Tat eine wesentlich kleinere Antinukleonenzahl, als die Theorie voraussagt. Die 
Abweichungen lassen sich qualitativ durch den Einflu8 der peripheren StoBe (siehe 
anten) und durch die Tatsache erklaren, daB ein Teil der Antinukleonen in dem- 
selben Kern vernichtet wird, in dem die Erzeugung stattfand. Das Verhaltnis 
ingx+):(nx-) sollte weniger als jenes von (np):(nx-) durch die Anwesenheit von 
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Kernmaterie beeinfluBt werden, d. h. es sollte auch bei Metalltargets einen Wert 
haben, der dem fiir pp-StéBe berechneten nahe kommt. Man findet experimentell 
bei einem K-Impuls von 2 GeV/c und einem Winkel von 18° (Laborsystem) (nix+)? 
(nx-) ~2 bis 3 (private Mitteilung von N. H. Lirmay, CERN), wahrend die 
Rechnung nach Transformation ins Laborsystem [25] folgende Zahlenergibt: Bei 
10° hat man 3.5, bis 20° findet man 3.85. Da beim Experiment auch p-n-St6Be bei- 
tragen, kann man verstehen, daB der gemessene Wert unter dem (fiir p-p) gerech- 
neten liegt. aes og 
Erst mit Wasserstofftargets wird man die 25 GeV-Rechnungen wirklich prifen 
konnen. 
DaB und welche Abweichungen dabei auftreten werden (die sich bei 6.2 GeV 
schon andeuteten) wird weiter unten diskutiert werden. 
4. Die pp-Vernichtung in Ruhe hat seit je Kopfzerbrechen gemacht, denn die 
beobachteten Pionzahlen waren um 5 verteilt. wahrend man mit der statistischen 


Exp. 300 6eV( Lab) ——> Yj 


Statistische 
erie a 


Experimental High energy Limit» L, 


1 
5; 
CUED (2) ab fe tab 
Tot, Athi eee CErye 
Bild 1. Bei 2.75 GeV und] 6.2) GeV liegen die ge- Bild 2. Bei 2.75 und 6.2 GeV hat man Ubereinstimmung 
rechneten Punkte praktisch auf den ge- zwischen Theorie und Experiment. Die Kurve 
messenen. Die Kurve ist yon Hand inter- ist durch diese Werte und zwischen den Grenzen 
poliert, so da’ sie durch alle gerechneten (je nach 2x) fiir 25 GeV gezogen. Der Wert fiir 
Punkte und die bei 25 GeV berechneten etwa 300 GeV (Emulsion) ist im Einklang mit den 
Grenzen geht (der Wert bei 25 GeV hingt gerechneten Werten. Der Wert bei 25 GeV 
noch von der Wahl 2x fiir K-Mesonen etwas (WINZELER, private Mitteilung) kann sich nur 
ab). Die Figur ist aus [78] enthnommen noch geringfiigig 4ndern, er ist gemessen worden, 


nachdem die Kurve publiziert wurde und hier 
nachtraglich in die aus [18] entnommene Ab- 
bildung eingezeichnet worden 


Theorie etwa (n,) ~ 3 berechnete. Verschiedene Autoren haben es mit Gewalt 
versucht: Man kann auf (n,) ~5 kommen, wenn man statt 2, = (42/8) (1/u)8 
ein mehr als 10mal gréferes 2, nimmt. Da aber bei der p-p-Wechselwirkung die 
Pionerzeugung gut wiedergegeben wird und dort ein zehnmal gréBeres Q, nicht 
tragbar ist, scheint diese Losung widerspruchsvoll. F. Cerutus [2/, 22] hat daher 
den von anderen Autoren (Literatur siehe [2/]) in anderem Zusammenhang 
vorgeschlagenen zx-z-Resonanzzustand n*(m ~4u,J =T =1) als instabiles 
Teilchen in die statistische Theorie der Paarvernichtung aufgenommen. Dann 
war es moglich, die richtige Pionzahl mit dem ,,natiirlichen‘‘ Wechselwirkungs- 
volumen (47/8) (1/)* zuerreichen. Auch die Spektren, der Anteil von K-Paaren [22] 
und die Haufigkeiten von 2: 4: 6-Spuren Ereignissen werden gut wiedergegeben. 
Das 2-Spuren-Spektrum hat am oberen Ende (~800 MeV kin. Energie des Pions) 
ein scharfes Maximum, herrithrend von dem direkt erzeugten xin p + p >7n*+ 2 
und ein weiteres Maximum von p + p->7 + x. Wahrend letzteres immer da ist, 
gibt es das erstere nur, wenn x* existiert. Seine Lage und Breite ist durch die 
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»,Masse‘‘ und Resonanzbreite (Lebensdauer) des ,,Teilchens 7* gegeben, Ex- 
perimentell liegen Anzeichen fiir ein solches Maximum vor, jedoch kann es sich 
um eine statistische Schwankung handeln; zur Zeit ist noch keine Entscheidung 
moglich. Das Pion-Isobar hat auch Winkelkorrelationen zur Folge. 

Da es hier nicht méglich ist, auf Einzelheiten einzugehen und viele berechenbare 
Gréfen experimentell kaum untersucht sind, sollen nur noch zwei Beispiele 
angefiihrt werden, die gewisse Ergebnisse der Rechnungen zusammenfassen [18] 
(p-p-StoBe) : 

In Bild 1 ist die GréBe [C — x+]/C als Funktion der Energie aufgetragen. C ist 
die Gesamtzahl geladener Teilchen, die das Reaktionszentrum. verlaBt (ein- 
schlieBlich der eingebrachten Ladung 2) und z+ die Zahl der geladenen Pionen. 
Demnach ist C — x+/CO das Verhiltnis der Zahl solcher geladenen Teilchen, die 
keine Pionen sind, zur Zahl aller geladenen Teilchen. 

In der nachsten Figur (Bild 2) findet sich C — 2 als Funktion der Energie. 
Neuere Vergleiche mit dem Experiment finden sich in [32] und [37]. In [37] han- 
delt es sich um die Produktion von Deuteronen, die befriedigend beschrieben wer- 
den kénnen. 


E. Diskussion 


Grenzen der Theorie und laufende Untersuchungen 


Der soeben gegebene Uberblick zeigt, daB die statistische Theorie in ihrer jetzigen 
Form manche experimentellen Fakten wiedergeben kann und das um so besser, 
je integraler die berechneten GroBen definiert sind. Absolute Teilchenzahlen (pro 
Ereignis) kommen ziemlich genau heraus, Spektren weniger gut. Offenbar sind 
auch die niedrigen (<2 GeV) und hohen (>25 GeV) Energien nicht giinstig: 
Bei den niedrigen spielen Details der Endzustande eine Rolle, die in der stati- 
stischen Theorie nicht behandelt werden (siehe aber eine Arbeit von B. R. Desay 


[26], in der Auswahlregeln fiir NN-Vernichtung mit der statistischen Theorie 
vereinigt werden) und bei den hdheren scheinen geometrische Effekte (StoB- 
parameter, Inelastizitat) wichtig zu werden. Diese letzten Effekte, mit denen wir 
uns jetzt etwas beschaftigen wollen, beruhen darauf, da bei wachsender Energie 
die Nukleonen im Schwerpunktsystem immer platter werden (Lorentz-Kon- 
traktion). Ihre relative Geschwindigkeit ist praktisch immer die Lichtgeschwindig- 
keit. Daher nimmt die Zeit, waihrend der sie miteinander im Kontakt sind, mit 
steigender Energie (d. h. mit zunehmender Abplattung) immer mehr ab. Die 
Wechselwirkung pflanzt sich aber senkrecht zur Bewegungsrichtung héchstens 
mit Lichtgeschwindigkeit durch das Nukleon fort. Wenn der Stof nicht zentral 
war, so kann es sein, da die Kontaktzeit zu kurz war, als dai die Wechselwirkung 
_sich iiber das ganze Volumen hatte ausdehnen kénnen. Dann trennen sich beide 
Nukleonen vor Erreichung eines Gleichgewichts und es entstehen zwei Zentren, 
die jedes fiir sich ein Gleichgewicht anstreben und dann einzeln durch die sta- 
tistische Theorie beschrieben werden sollten [27]. Das wiirde zwei Effekte erklaren : 


— Erstens wird wegen der zu kurzen Kontaktzeit die Geschwindigkeit nur teil- 
weise abgebremst und in Anregungsenergie verwandelt. Jedes Zentrum behalt 
noch einen unter Umstinden betriachtlichen Anteil kinetischer Energie und 
nur den Rest kann als Anregungsenergie zur Produktion von Teilchen benutzt 
werden. Die Inelastizitat wird also bei fester Energie mit wachsendem StoB- 
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parameter und bei festem Stofparameter mit wachsender Energie kleiner 
werden. Damit wird die Baryon-Paarerzeugung gedriickt. 

— Zweitens bewegen sich die beiden Zentren. Wenn jedes in seinem Ruhsystem 
isotrop emittiert, so ergibt sich im gemeinsamen Schwerpunktsystem eine 
Vorwiairts-Riickwarts-Anisotropie. 


StéBe bei denen — bei gegebener Energie — der Sto&parameter so klein bleibt, — 
daB die Wechselwirkung sich waihrend der Kontaktzeit iiber beide Nukleonen 
ausdehnen kann, diirfen mit der vorliegenden 
Theorie als ,,Zentrale StdBe‘‘ behandelt werden. 
t Man sieht, dafB die Definition von ,,Zentral“ energie- 
6 abhingig ist. Nach obigen Uberlegungen liest man 

| aus Bild 3 ab: 
Der Nukleon-Durchmesser ist d ~ 2/yu mit 
y = (1 — fp)" = E/2, wenn £ die totale Energie 
im Schwerpunktsystem ist (MZ = 1!). Die Kontakt- 
zeit ist also t ~2/yu, und wenn b >t ist, kann die 
Wechselwirkung nicht beide Nukleonen ganz er- 
fassen. Man hat also die Bedingung fiir ,,zentralen 


ase StoB 
ye 4 
Bild 3. Zur Definition ,,zentraler*‘ b = bo SS SS Se (42) 
StoBe yu we 
Da der Nukleonendurchmesser 2/4 das natiirliche Maf8 ist, finden wir fiir 
Ozentral == 6? ? < (=) (43) 
Ogeometr. (2 / 1) E 
das ergibt folgende Werte: 
Ftep,kin (GeV]| 275 | 62 | 25 
BE ERY ere OEM 8 ee 
07,/z, | « 0.4 | 0.23 | 0.07 


Diese Uberlegungen zeigen, da8 man versuchen sollte, bei gegebener Energie nur 
die ,,zentralen‘‘ St6Be mit der vorliegenden Theorie zu beschreiben. Dann muB 
man fiir die gréBeren StoBparameter ein Zwei-Zentren-Modell [27] benutzen. Die 
Anregungsenergie 4H wiirde man etwa proportional dem Volumen setzen, iiber 
das die Wechselwirkung sich wahrend der Kontaktzeit ausbreiten kann. Dann 
ware 


AV 


AH ~—,--(H—2)=K(b, BE): (E—2). (44) 


(Nukleon-Nukleon-StoB) 


Der Rest, H — 2 — ALF, ist die kinetische Energie der beiden Zentren, womit 
(unter der Annahme, dafB keine nennenswerte Richtungsanderung stattfindet) die 
Transformation vom Ruhsystem der Zentren auf das Schwerpunktsystem voll- 
standig und die Winkelverteilung fast vollstandig festgelegt sind. AE|(E — 2) 
ist im wesentlichen die Inelastizitiit. 
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Kin einfaches Modell kann qualitativ zeigen, was dabei zu erwarten ist. Nehmen 

wir an, die Nukleonen seien homogene Kugeln!) mit dem Radius r — 1 [u. Wir 

bezeichnen mit AV das gesamte (eng- und weit-) schraffierte Gebiet in Bild 4 

und bérechnen AV unter der vereinfachenden Annahme, daB dieses Gebiet eine 

eee der Hohe a = 2r—b4b, sei. Aus elementaren Formein folgt 
ann 


AV b— b,]2 b= 6 
ee [1 a] [t+ | fir 6<6b<2r (45) 
1 fur VR) = 6; 
4r 
6, = ak r = — = ,,Nukleonenradius‘. 
lt 


AE = (E — 2): K(b, E) ist die zur Teilchenerzeugung angebotene Energie. Wir 
nennen K (b, #) die ,,angebotene Inelastizitat‘‘. Fiir hohere Energien wird sie der 
experimentellen Inelastizitat nahe kommen. 


Nachdem K(b, Z) durch ein Modell festge- K 
legt ist, kann man bei gegebener Energie H die b 
Wabrscheinlichkeitsverteilung ; i 
W(K, B)dK = logis 
r E FR Ss fr 
— 22 b(K, E) oe \ aK (46) Lele 


durch Umkehrung von (45) berechnen. Da 
mit K(b, H) die kinetische und Anregungs- 


Bild 4. Zur Inelastizitaét: Definition von 4 V/V. 


energie der Zentren festgelegt ist, kann man Das*Snavelinatierte Gebiek Yat daw Ace 
auf jedes Zentrum (etwa unter Annahme Uberlappung wihrend der Kollision; 

2 = Taxis das weitschraffierte ist das Gebiet, in 
der Gleichverteilung der Anregungsenergie: ag be ae ee eee Se 


AE/2 je Zentrum) die statistische Theorie Kontaktzeit 27/y eindringt 
anwenden und auf das gemeinsame Schwer- 

punktsystem zuriicktransformieren. Das ergibt fiir jedes 6 eine Impuls- und Winkel- 
verteilung, Multiplizitat usw. Zuletzt hat man iiber alle b (bzw. nach (46) tiber alle 
K) zu integrieren, kann aber auch mit (46) Wahrscheinlichkeitsverteilungen an- 
geben. Mit der in dieser Arbeit beschriebenen Technik scheint das fiir héhere 
Energien aussichtslos, da die erforderlichen Rechnungen zu lang werden. Bei 
6.2 GeV ist eine derartige Rechnung mit Erfolg durchgefiihrt worden [20]; wir 
kommen darauf zuriick. 

Wir wollen aus (45) nur eine Folgerung ziehen, nimlich die (iiber alle Stof- 


parameter b) gemittelte ,,angebotene“ Inelastizitat K (EB) berechnen, Offenbar ist 


2r 


K (£) = aaa | 22? db- K(b, BE). (47) 
0 


1) Natiirlich sollte man, wenn dieses Modell zu quantitativen Rechnungen benutzt wird, die 
Nukleonenstruktur beriicksichtigen. 
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Man erhalt mit (45) 


il 4r 


Agr? 


2r 
K (BE) =zal bo? + [ 2b ab- KO, B) | bo= mae (48) 
bo 


Man sieht, da® K(E) eine Funktion ist, die von H = Epin = 2 an stetig abfallt. 
Sie erreicht einen konstanten Wert, wenn EH -> oo, also b,—> 0 geht, da dann nur 
das Integral iibrig bleibt. Man erhalt aus (45) 


ac 1\2 4 5 | 
E(z) = 5+ oltte[(r-Z) +l) B>2. (49) 


Fir kleine Energien sollte man K(E) nicht mit der experimentellen Inelastizitit 


vergleichen, bei sehr groBen Energien erhilt man K (HZ) = 0.15 + 1/H, was der 
in kosmischer Strahlung gefundenen Verteilung nicht widerspricht (v. LINDERN, 
private Mitteilung). 

Wir wollen die Sache hier nicht weiter verfolgen. Es sollte nur angedeutet werden, 
da8 an dem naiven Bild doch etwas Richtiges sein méchte. 

Die soeben gegebenen Uberlegungen zur Inelastizitét als Funktion von Energie 
und StoBparameter kléren auch einen weiteren Punkt auf, namlich die Dis- 


krepanz zwischen N-N-Sto8en und N-N-Vernichtung in Ruhe. 


Wir hatten erwahnt, daB mit Benutzung eines z-z-Resonanzzustandes die N-N- 
Vernichtung in Ruhe in Ubereinstimmung mit dem Experiment beschrieben 
wird [21, 22]. 

Konsequenterweise sollte dieses x*-,,Teilchen’’ auch bei N-N-Wechselwirkung 
eingefiihrt werden. Dann werden, wie eine Rechnung bei 6.2 GeV zeigte, die 
Multiplizitaten hoher als die experimentellen. Eine Beriicksichtigung der nicht- 
zentralen StoBe (nach der Tabelle hinter (43) also ~76°) mit ihren kleineren 
Inelastizitaten wiirde das wieder ausgleichen und zudem eine Winkelverteilung 
liefern. Fiihrt man das (45) zugrunde liegende Modell quantitativ durch [20], so 
ergibt sich, dai fiir die nichtzentralen St6Be die Anregungsenergie gerade aus- 
reicht, um die sich bewegenden Zentren in den ersten N-z-Resonanzzustand zu 
versetzen, also entweder N*-+ N* oder N + N* oder N + N zu erzeugen. Die 
Autoren vernachlassigen alle anderen Anregungszustinde in diesem Energie- 
gebiet (der Erfolg des LinDENBAUM-STERNHEIMER-Modells [6] bei <1 GeV 
Primarenergie rechtfertigt das) und fiihren ferner das Verhaltnis ozentrai/Ogeom als 
Parameter ein (der die GréBenordnung 0,23 haben sollte). Multiplizitaét, Haufigkeit 
der Ereignisse mit 2, 4,6 Spuren, mittlere Energien und Winkelverteilung werden 
berechnet. Die experimentellen Ergebnisse [28, 29] stimmen untereinander nicht 
vollig tiberein. Daher ist in Einzelheiten noch kein Vergleich méglich, doch kann 
man die Zahl der geladenen erzeugten Teilchen und die Inelastizitat mit ozentr |Ozcomt 
~ 0.5 richtig wiedergeben. Damit ist der einzige freie Parameter bestimmt. Das 
Verhaltnis der Ereignisse mit 2, 4,6 Spuren kommt dann richtig heraus, ebenso 
die mittleren kinetischen Energien. Dabei scheint die Anwesenheit eines Pion- 
Isobars x*(J = T’ = 1; m ~ 4) wesentlich zu sein. 

Diese explizite Rechnung, die leider mit den z. Zt. vorhandenen Methoden kaum 
fiir héhere Energie durchfiihrbar ist, zeigt, daB hier wenigstens eine Moglichkeit 
liegt, weiterzugehen. Hs ist eine Frage der technischen Durchfiihrung. Die hier 
beschriebene Methode bleibt ganz im Rahmen der statistischen Theorie. Eine 
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wesentlich detailliertere Behandlung, aber auch noch im Rahmen der statistischen 
Theorie findet sich in [23, 24]. 

SchlieBlich kann man versuchen, die peripheren StéBe ganz aus der statistischen 
Theorie zu lésen und mit den iiblichen Methoden der Pionphysik zu behandeln 
[30]. Dann ist die Zusammenfiigung der beiden Methoden eine noch ungeloste 
Aufgabe. 

Nach dem Gesagten scheint es merkwiirdig, daB® bei 25 GeV die durchgefiihrten 
Rechnungen [18] nicht véllig falsche Ergebnisse liefern. Entweder ist unsere 
Schaétzung Von ozentrai/Ogeom, ZU pessimistisch — worauf das Ergebnis der Unter- 
suchung bei 6.2 GeV [20] hinweist1) — oder die Behandlung aller St6Be als zentral 
(wie in [18] geschehen) itberschaétzt die Multiplizitiiten in solcher Weise, daB 
damit die Vernachlassigung des Pion-Isobars x* und hoherer Anregungszustande 
der Nukleonen kompensiert wird. Vermutlich wirken sogar beide Ursachen ge- 
meinsam. Hier bleiben, teilweise wegen der schwierigen experimentellen Situation, 
viele Fragen offen. 

Der vorliegende Bericht beschreibt im Wesentlichen die bei CERN in den letzten 
beiden Jahren gewonnenen Ergebnisse. Wir wollen — ohne ins einzelne zu gehen, 
was den Rahmen dieser Arbeit sprengen wiirde — wenigstens erwahnen, daf 
seit der Entwicklung von Beschleunigern, die in die Gegend von 10 GeV und mehr 
vordringen, an allen Orten, wo solche Maschinen stehen, an dhnlichen Fragen ge- 
arbeitet wird, vor allem in Dubna (UdSSR). 

Meinen Kollegen bei CERN, besonders den Herren F. Cerunus und J. von Brnr 
bin ich fiir die enge Zusammenarbeit dankbar. 
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Inelastic Interactions between High Energy Particles") 


V.S. BARASHENKOV 
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Up to now for the description of inelastic collisions between high energy particles 
the Fermi theory was used. This theory is based upon a statistical equilibrium 
between the particles generated in the collision. A great number of papers are 
devoted to the development of this theory (the bibliography is given, e. g., in [7)). 
However, the accuracy of the results obtained was not high as the numerical 
calculations are very complicated. Only recently mathematical difficulties have 
been overcome, and new effective methods of calculations have been developed. 
The application of electronic computers allows to make the calculations quickly 
and effectively up to the energies of order of some dozen Bev. On the other hand, 
a great deal of experimental data obtained by means of the USA and USSR 
accelerators made it possible in some cases to compare the experiment with the 
theory. It has been thereby established that in the framework of the Fermi 
statistical theory itis impossible to account for many experimental facts already 
known. For their understanding it was found necessary to take into account the 
intrinsic structure of elementary particles. 

Let us consider these problems in more detail. 


1. Comparison of Experimental Data with Calculations by Fermi 
Statistical Theory 


It has been shown in the papers of many authors that the statistical theory of 
multiple production accounts for the average number of generated charged and 
neutral particles fairly well. The theoretical distributions by the number of prongs 
of the stars which are due to the collisions of high-velocity pions and nucleons 
with nucleons and in the antinucleon annihilation are found to be very close to 
the experimental ones. In Fig. 1, as an example such a distribution is given for 
the interaction of 6.8 Bev pions with protons and neutrons. 
Close are also the experimental and theoretical distributions of generated particles 
by the momenta. In Figs. 2 and 3 are compared the experimental and theoretical 
distributions for the (x-N) interactions at H = 6.8 Bev. (Experimental data are 


taken from [4]). 


1) A review at the All-Union Inter-University Conference on the Field Theory and the Theory 
of Elementary particles; Uzhgorod, 1960. 
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However, in all cases the experimental spectrum of nucleons turns out to be 
harder than the theoretical one. For instance, for (PP) collisions at 9 Bev the 
experimental value of the mean momentum of protons in the c. m. s. is found to be 
(1.24 -- 0.25) Bev whereas the theoretical one is 0.8 Bev. As for the mean 
experimental value of the energy necessary for the production of new particles per 
collision, it turns out to be smaller than the calculated one. 

The statistical theory of Fermi contradicts the experiment sharply when the 
angular distributions of generated particles are compared. In accordance with the 
main idea about the statistical equilibrium, which is the basis of this theory, the 


0 7 2 3 4 5 6 7 én 


Fig. 1. es ratio of the number of n-prong stars experimentally observed to that calculated by the statistical 
heory. 
An odd number of prongs —(z n) collision 
An even number of prongs —(m p) collision 
Energy H = 6.8 Bev. 
The errors in the experimental data are indicated 


angular distribution of the generated particles must be symmetrical in the c. m. s. 
with respect to the angle 0 = 2/2. The account of the law of momentum con- 
servation leads to the anisotropic angular distribution of the generated particles. 
For example, in case of two-particle reaction this distribution is as follows [22]1): 


oF 2 li (M2 — M2) fea ong 
dQ V¥ Ei EORey 


EB? R? 
x a rues — 14) 8(Cos 0), 


1) This expression passes into a well-known BLook’s formula [25], if it is i 
anda Tee mula [25], if it is integrated over dQ 
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where : nt 
Y ,_ 20a te) es 
+. Ey Es i 
Gin Oe, 
2M h : 
E, is the total energy in thec.m.s., R = ee ; M is the nucleon mass, 
Cc 


0 . . . 
pis the mass of a x-meson; M, and M, are the masses of the colliding particles (in 
the general case M;-+ M and M,=- m); L is the momentum of the colliding 


particless; 6 isthe angle between the vector of the angular momentum L and the 


incident particle velocity; 
6(z) =1, when ~>0, and 6(z)=0, when x<0. 


Obviously, 
dW(o) dW(x—8) 
Co. dQ : 


However, the angular asymmetry in the Fermi theory can be obtained in no way. 
The experimental angular distributions are found to be very asymmetrycal. For 
instance, in case of (xp) collisions at H = 6.8 Bev the ratio of the number of 
protons emitted in the front half-sphere n, to that of protons emitted in the back 
half-sphere n, (with respect to the velocity vector of the primary ~~-meson) is [4]: 


Ay =Mo/% <KO1L01 << 1. 


An analogous ratio charged pions in this case is equal to 1.57 + 0.1 [4]. 

These experimental result point out that the collisions are possible when there is 
no statistical equlibrium between all the generated particles. In these cases a 
more detailed consideration of the collisions is necessary. 


A Model of Central and Peripheral Collisions 


An analysis of elastic interactions of different particles with nucleons makes one 
to think that the matter in the cloud of virtual particles surrounding the centre 
of the real physical nucleon is distributed very inhomogeneously (non-uniformly) ; 
in the central regions its density is much greater than in the periphery. The 
particles in the central regions of a virtual cloud are considereably stronger bound 
with each other, than with the particles from the peripheral layers of the cloud [4]. 
Such a model of a nucleon suggests three types of nucleon-nucleon collisions: 

I. The collision between the central nucleon regions, when in a small space volume 
the major part of the energy of the colliding nucleons is released, and an intensive 
production of new particles is taking place. 

II. The collision of the periphery of one nucleon with the central part another 
nucleon. So to speak, in this case an “‘stripping” of the peripheral mesons in a 
nucleon occurs. 

IIT. The collision between the peripheral regions of the nucleons. Since the density 
of matter on the periphery is much smaller than in the central regions of a nucleon, 
such collisions will be considerably rarer than two other types of collisions. 

Such a classification of (NN)-collisions was first discussed by D. I. BLoKHINTSEYV 
[6]. This classification, of course, is applicable only for the collisions of high 
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velocity particles with the wave length 4, much smaller than the dimensions of 
these particles. Only in this case one can speak about the collision parameter. 
Practically this holds already for particles with E ~ 1 Bev. 
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ig. 6. The angular asymmetry of protons Ap = itp | Np in the c. m.s. of the colliding 7 -meson and proton. The 
region of possible values Ay is shaded. The horizontal shading is the variant with the isobar, the vertical 
shading — the variant without the isobar 


\t present there are also theoretical and experimental indications that the 
iteraction between x-mesons is strong. This means that around the center of the 
sal physical meson there must exist a cloud of virtual ~-mesons. Since the 
-mesons strongly interact with nucleons, then inside a physical meson there must 
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also be a cloud of virtual nucleon-antinucleon pairs. It is possible that there 
exists a cloud of virtual K-particles etc. In other words, a real physical x-meson 
must have dimensions and intrinsic structure. From the standpoint of this 
model (xN) collisions may also be divided into three types like (NN) collisions 
what was done above. 

An analogous classification is valid also for the collisions of high velocity particles 
of other kinds. 

The account of peripheral collisions leads to an anisotropy of angular distributions ~ 
in the c. m. s. of the colliding particles (see Fig. 4). If the angular distribution of 
particles, generated in the collision of the particles 6 and c, is close to an isotropic | 
one in the c. m. s. (b; c), then in passing to the c. m. s. (a; ¢) this angular distribu- 
tion will be sharply anisotropic due to the compression of angles in the coordinate 
transformation. 

However, in some cases, e. g., for (pN)-collisions, the angular distribution in the 
c.m. 8s. of the colliding particles becomes anisotropic, the asymmetry, however, is 
not obtained [2]. In [7, 2] it has been shown that the angular asymmetry in this 
case may be obtained only by assuming that the nucleon which emits the peri- 
pheral meson passes into an excited isobaric state, which, in its turn, decays in- 
dependently (see Fig. 5). An analogous model was suggested by I. E. Tama [7]. 
The calculations have show that all the published experimental data on (NN) and 
(x N)-collisions may be accounted for if the cross section for the peripheral collisions 
amounts not less than 20°,—30% of the total cross sections for all inelastic 
processes [2, 8]. However, the analysis of more accurate experimental data 
recently obtained at Dubna on (z-N) interactions at H = 6.8 Bev (see [4]) 
shows that the cross section for the peripheral collisions o, must be, at any rate, 
greater than the half of all the inelastic (xN) interactions. This is seen from Fig. 6, 
where the theoretical angular asymmetry of protons A, is plotted against the 
parameter € = o,/ojn. (Let us remind, that the experimental value is 


Ay 0 ee Ot 


When = 0.5 the experimental data on (NN) collisions are in better agreement 
with the theory, than with § = 0.2 = 0.3. 


3, Energy Spectrum of Peripheral Mesons in a Quick Nucleon 


At present there is a series of preprints whose authors treat different problems of 
the theory of peripheral collisions. However, this theory is still very incomplete. 
In particular, for the calculations the energy spectrum of the peripheral mesons in 
a nucleon is usually used a formula which has been obtained by Herter and 
PENG [9] many years ago. This formula was obtained under the assumption that 
the meson field in a nucleon is described by the classical expression 


| ite d ere 
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(h or 1; w is the mass of a z-meson). However, this expression, as was shown 
in [4] is a very bad approximation in the region r < 1/y. 

It is possible to get a more accurate expression for the energy spectrum of the 
peripheral mesons [24]. 
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In the general case this spectrum is determined from the comparison of two 
expressions for the total flux of the meson field energy of the nucleon in a layer 
with the radius 0 (see Fig. 7) 


Cc 


I(o) =f eq(gse)de 
and. > 


I(9) =2x%@ [ S(g;t) dt. 
+c 


Here q(o;¢) is the energy spectrum of mesons at a distance o from the center 
of a quick nucleon; S(g;¢) is the 2-th component of the Poyting-Umov vector 
at the point « —0 at the moment of time ¢ (see Fig. 7). It is easy to show that 


bowels (l., — Ty (1 6) Td, 


where 7';; is the mean value of the component of the meson field energy-momen- 
tum. tensor in the nucleon coordinate system; y? = (1 — f?)~*/:; B is the velocity 
of a nucleon in the Lab.system (parallel to the X-axis). 

As we are considering the periphery of a nucleon, then in the first approximation 
the nucleon recoil may be neglected. In this case for calculating 7';; it is possible 
to apply the CoEw-Low theory [10]. In the one-meson approximation 
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where J, and J, are the well-known Bessel functions; v(w*) is the form-factor of 
the meson field source; wz = k? + e2/y2f2 + w?, w is the mass of a pion. The 
total energy spectrum of the peripheral mesons in the nucleon shell with the radius 
eo fis 


g(e) = 2x f 9 q(o; €) dee 
R 


Now a numerical calculation of this spectrum is being made by means of an 
electronic computer. 
4. Resonance Interaction of ™-Mesons 


Another problem which attracts a great deal of attention is the account of the 
possible resonance (x7) interaction and its influence on the calculations according 
to the statistical theory. From the theoretical point of view the idea of such an 


3* 
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interaction seems very attractive and plausible. However, all the known experi- 
mental data on the multiple production of particles and their energy spectra in 
(xN) and (NN) collisions are possible to be explained now without taking into 
account there resonance (x7) interaction. In the framework of the Fermi theory 
the account of such an interaction makes the calculations much more difficult and 
does not give a better agreement with the modern experiments. For example it can 
be seen from Table I and II, where are given the experimental and theoretical 
data for (x-P) interactions at H = 5 Bev and (PP) interactions at 6.2 Bev (ex- 
perimental data from [11, 12]). The calculations have shon that an analogous 
situation occurs in (x-N) collisions at H = 6.8 Bev and in (NN) collisions at 
Hi = 9 Bev?). 

In paper [13] it is stated that such conclusions are obtained because the statistical 
theory is crude. But this just means that within the accuracy of the present 
theory and experiment there is no necessity of taking into account the resonance 
(x7) interaction in the analysis of (xN) and (NN) collisions. This concerns 
especially the energy region near 1 Bev, where the statistical theory is applicable 
with reserve. 

On the other hand, an account of the resonance interaction of z-mesons con- 
siderably improves the agreement of the theory with experiment in case of a slow 
antinucleon annihilation. This is seen from Table III, where are given the results 


Table I. 
(x-p) interaction; H = 5 BeV 


Theory 
Chavakierisie Without With account | With account 
Bflatcolliston A keounbio? of (x7) inter- | of (x77) inter- | Experiment 
Gen antes action; action; 
& Aeiar Dyson’s TAKEDAS’ 
variant”) variant?) 
The number of stars 
with 2-prongs [in %%] 61,5 50 44 60,8 + 7,5 
ye % 9 37 44 51,5 36,4 + 5,8 
2? 6 ” 2 ” 1,5 6 4,5 2,8 + 1,6 
The average number of all 
charged particles 2,8 3,1 3,2 2,84 + 0,48 
The average number of all 
charged and neutral 
particles 4,2 ANT 5,1 4,3 


1) The results of the calculations are close to the data obtained by HaGEDoRN and CERuLus 
using the Monte-Carlo method [18]. I am grateful to Drs. Hagedorn an Cerulus, who sent ee 
results of their calculations. 

*) In Dyson’s variant the spin of a z-meson isobar is S = 0; its isotopic spin is 7 = 0. In 
TAKEDA’Ss variant S = 0; 7’ = 1. In both cases the mass of a x-meson isobar was easel to 
be equal to three z-meson masses: y* = 34; the change of u* from 3u up to 4 and the spin 
from S=0 upto S=1 does not strongly affect the results of the calculations. : 
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of the calculations of the stars due to the annihilation of slow antinucleons, as well 
as are listed corresponding experimental data from [15]. 

It is also possible to bring into agreement the theory and experiment without 
resorting to the hypothesis about the resonance (x7) interaction. However, it 


appears necessary to make some special assumptions, the validity of which is not 
quite obvious (see, e. g. [14, 15]). 


Table II. 
(pp) interaction; # — 6.2 Bev 


a 


Theory 
Characteristic of a collision ieee a With account Experiment 
| (m7) inter- of (rm) ae 
action ett 
a 
The number of stars 
with 2 prongs [in 9%] 43 26 32521221 
2 53 64 53,8 -+ 10,7 
poh ae 5 os 4 10 11,8 -- 2,6 
pO s ees 28s 0,1 0,2 11,8 + 11,3 
The mean number of neutrons 0,8 0,8 = 
2 >» >, protons ie 12 1,6 -—- 0,6 
- in Fe ;, 7-mesons 1,0 1,2 _ 
a zs - ;, ™=+-mesons 2,0 2,5 1,9 + 0,3 
ws z 55 ,, all charged 
particles 3o2 ahi 3,5 + 0,3 
The mean momentum of neutrons 
(c. m. s. Bev/c) 0,78 0,71 = 
[The mean momentum of protons 
(c. m. 8. Bev/c) 0,78 0,71 0,77 + 0,01 
[The mean momentum of 7°-mesons 
(c. m. s. Bev/c) 0,51 0,41 _ 
[The mean momentum of 7+-mesons 
(c. m. s. Bev/c) 0,51 0,41 0,32 
[The mean momentum of all charged 
particles (c. m. s. Bev/c) 0,61 0,52 0,53 


) The mass of a z-meson isobar u* was assumed to be equal to four masses of a 7-meson; 
= 4u; 5S =i; T= 1. 


Thus, in the framework of the statistical theory a succesive consideration of the 
esonance (x7) interaction is rather difficult. However, an account of the peri- 
yheral collisions changes the situation essentially. In peripheral collisions a some- 
vhat less number of particles is produced than in central collisions (see Table IT 
rom [8]). We hope this will compensate the increase of the multiplicity due to the 
ccount of the resonance (x7) interaction. This will be clear when detailed 


umerical calculations are completed. 
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Table III. 


(NN) 
(pp) 


interaction 


Theory 


With account 


With account 


nie 


Characteristic Without 3 : 
of a collision account of of (mz)-inter- | of ngs Experiment 
Gale action action Ba 
action aaa tk 
S=1; 7=1|/S=1;7=1 


The number of stars 
with 2 prongs in [%%] WOH 27 45 42 +7,3 
rey Agar ans ames 20 65 52 52 + 8,1 
yb GON es ls ees 0,1 8 3 6,3 + 2,8 
The average number of 
7°-mesons de 1,7 3s 1,60 + 0,50 
The average number of 
™+-mesons AD! 3,6 1,6 3,06 + 0,16 
The average number of 
all charged and neutral 
particles 3,4 5,3 4,7 4,7 +0,1 
The mean momentum of was not 
m™-mesons (in Bev/c) calculated 0,33 0,35 0,342 + 0,05 


5. Multiple Production of Strange Particles 


At the present time still very incomplete experimental information is available 
on interactions between strange particles at high energies, and, in particular, on 
the multiple production of these particles. All what is known is that the cross 
section for the generation of strange particles is only some percent of the total 
cross section for all the inelastic processes. 

To account for this experimental fact, in paper [16] a statistical model was 


suggested, where the multiple production of K- and K-particles takes place in a 
smaller space volume than of ail other particles. Irrespective of its physical inter- 
pretation such an approach has become now generally accepted (see e. g., [17, 18]). 
At the same time, as more and more experimental data become available, such a 
simplified approach is no longer satisfactory. 

In particular, the calculations on the production of strange particles in the 
annihilation of slow antinucleons lead to too low values of the cross sections 
compared with the experimental ones (see, e. g. [17] or [13]). The experiments 
with (x~N) interactions performed at Dubna at an energy E = 6.8 Bev also 
indicate that the number of the generated K-particles is found to be larger than 


that calculated theoretically under the assumption that K- and K-mesons are 
produced in the same volume with the radius rg = i/Mxc = 0.4: 10-13 cm, less 
than that for the production of particles of all other kinds. However, the total cross 
section for the production of all strange particles is found to be, very likely, 
somewhat greater than that calculated theoretically at rg—= 0.4 + 10-13 em. 

These results indicate, that by this time a more detailed consideration of the 
mechanism of strange particle multiple production is necessary. The introduction 
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ofa smaller space volume is only a means of the roughly phenomenological account 
of different strength of the interaction between different particles, This approach 
is no longer satisfactory. 

a further progress in the theory more detailed experimental data are necessary 
also. 


6. Interaction of Particles at Very High Energies 


Cosmic ray experiments allow to study the interactions between particles up to 
gigantic energies H~ 10° Bev. These experiments by their nature may be divided 
at present into two groups: 

1. A study of wide atmospheric showers due to the primary particles from the 
cosmic space. 

2. A study of the tracks of the cosmic particles interaction in the emulsion stacks. 
In both these cases to get an idea on the interaction between elementary particles 
is a very complicated and often an ambiguous problem: in the first case the 
character of the primary act of the interaction is extremely obscured by the 
interactions of successive generations of the produced particles; in the second 
case the interaction of a cosmic particle takes place with a nucleus, and the features 
of elementary (xN) and (NN) interactions also turn out to be spread. Besides, 
the statistics of the recorded events is usually poor. 

Nevertheless, a detailed analysis of numerous experiments still allows to obtain 
some conclusions: 

1. The cross sections for (<N) and (NN) interactions remain constant, or, at any 
rate, do not increase up to EH ~ 10° Bev. 

2. The total number of secondary particles increases with the energy of the primary 
cosmic particle # quicker than 1/3. At the same time an overwhelming majority 
of these particles are z-mesons. 

3. Among the generated particles there is observed, as a rule, a particle with a 
very great momentum. Such a particle possesses a large part of the energy of a 
primary particle. (According to the data of different authors — from 30% up to 
70%). In the center-of-mass-system of the colliding particles the transverse 
momenta of secondary particles are found to be ~0.5 Bev/c irrespective of the 
emission angle of these particles. 

4, The angular distribution of the generated particles in the c.m.s. is sharply 
anisotropic. At the same time in many analyzed cases the angular distribution is 


found to be asymmetrical: A =nln = ae 

To account for the experimental data a hydrodynamical theory of multiple pro- 
duction is usually made use of. Recently, however, other models were suggested 
(for instance, a model of a “‘fiery ball” ete.). An analysis of the experimental and 
theoretical results is given in review [20]. Let us point aut to two problems which 
are of interest in connection with the interpretation of the experiments with the 
accelerators. 

At high energies of the colliding particles an account of the peripheral inter- 
actions is of still more importance than at those obtained by accelerators. The 
asymmetry experimentally observed in the stars A =- i” and especially the role 
of one of the secondary particles qualitatively confirm this conclusion. From this 
standpoint, the results obtained by means of the hydrodynamical theory and in the 
first turn the calculated angular distributions of secondary particles require a 


critical treatment. 
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The second question concerns the interaction of quick particles with nuclei. 

At energies up to 10 Bev the experimental data on nucleon and meson inter- 
action with nuclei are in good agreement with the mechanism of the intranuclear 
cascade. An agreement between the experiment and the theory holds both for ~ 
light and for heavy nuclei. The calculations have shown [21] that to account for 
the multiplicity and energy distributions of quick and cascade particles it is 
sufficient to take into account two collisions in light nuclei, and three collisions 
in heavy nuclei. The angular distributions of quick particles are found to be close 
to the experimental ones. To account for the angular distributions of the cascade 
particles it appeared necessary to take into account one more collision. 

At energies HS 10 Bev both the incident primary particle and the nucleus (in 
the c. m.s.) are strongly compressed. Therefore, the time of interaction with a 
nucleus is very short and the perturbation wave has no time to expand in a 
transverse direction further than by h/wc. The interaction under these conditions 
occurs only with the tube of the nuclear matter cut by the primary particle. 
However, even at very high energies, the mechanism of this interaction, in 
principle, may be different: either like at lower energies it will be a cascade inter- 
action with separate nucleons, or, as a result of the interaction a unified excited 
“compound-system” is formed, which further decays into separate particles. 
The latter mechanism has been recently investigated in detail in [23]. It has been 
shown that for the most part the theory agrees qualitatively with experiment. 
However, a quantitative difference often amounts to 30%—40% and more. 
Within experimental error, the results of the measurements may be also brought 
into agreement with the theory of the intranuclear cascade. 

Thus, the mechanism of high energy particle interaction with atomic nuclei 
remains at present still very obscure. First of all, more careful experimental 
measurements are needed. 
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Small Angle Elastic Scattering of High Energy Particles’) 


V.S. BARASHENKOV 
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In theoretical investigations it was pointed out long ago that elastic scattering of 
high energy mesons and nucleons is an effective means for the study of the 
intrinsic structure of elementary particles and atomic nuclei (see, e. g. [1 —4], 
where a detailed bibliography is given). 

In particular, the small angle measurements allow to get the information on the 
peripheral interactions of particles. 

Some problems associated with the interpretation of the experiment on elastic 
scattering of quick particles will be treated further. 


1. Real Part of Elastic Scattering Amplitude at High Energy 


The experimental data obtained recently at Dubna on elastic (pp)-scattering at 
8.5 Bev have shown that the model of a purely absorbing nucleon so far widely 
used is too rough and as it seems, cannot account for all the experimental results. 
The magnitude of elastic scattering cross section at zero angle 6 =0 from [5] 
greatly exceeds the value obtained by using an optical model (assuming that the 
total cross-section for (pp)-interaction at 8.5 Bev o ~ 30- 10-2? cm? [5]). The 
experimental data may be formally accounted for by assuming that the real part 
of the elastic scattering amplitude D(H) at 0 =0 is different from zero and by 
the module is equal approximately to 3-10-% cm (in c. m. system). An unex- 
pected and surprising result is that D(H) takes on so great values at HS 1 Bev. 
In this connection a theoretical investigation of the behaviour of D(H) at high 
energies is of great interest. 

In [6] it has been shown that at H — co D(£) increases at any rate more slowly 
than # In #. However, it is possible to obtain more definite conclusions about the 
behaviour of D(£) at high energies. 

To this end we put the elastic scattering cross-section og] as: 


co co 


oo =H 12D) (QU + 1) [1 — fy eR = ow A (20-4 1) (1 — By)? + 
1 


1=0 =0 


co 


+ 4a #2 S (21+ 1) fi Sin? §) = 3+ Dp - 


=0 


1) A review at the All-Union Inter-University Conference on the Field Theory and the Theory 
of Elementary Particles, Uzhgorod, 1960. 
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Here the first term is wholly due to inelastic processes (6; = 1) and is a cross- 
section of a purely diffraction scattering. 
The second term depends upon the non-diffraction scattering. 


When E>1 


D2 © 821? [ B,157d1 = const. B (E) - 62 (B) , 
0 


where 9 = // =1/E is the collision parameter, whereas §(Z) and 6(£) are values 
of the functions in the range where the nuclear forces act. 

Since at high energies all the scattering becomes a purely diffraction one, then 
6(£) — const. and 6(#) — 0. It follows from here, that 


D(E) =24 3" (21 + 1) B, Sin26, ~ const. H - 6 (BE) (1) 
t1=0 


increases at any rate, more slowly than Z. 

This result means that in the dispersion relation for D(#) all the term, whose 
degree of increasing greater or equal to #, must vanish. From here it follows that 
in the particular case of elastic scattering of pions on nucleons 


po — pe SP. Bf AE fo,(B) — 0B) =O, @) 


“ 


where D*. = D..(u); uw is the pion mass; k = {E? — uw; ox are the total cross- 
sections for x*-meson interaction with protons [7]. 

Making use of the experimental values of the cross-sections o from relation (2) 
it is possible to determine the value of the meson-nucleon coupling constant /? by 
means of the statistical analysis. The calculations have shown that?) 


2 — 0.08 + 0.005. 


Tn case of elastic nucleon-nucleon scattering the relation analogous to (2) is of the 
form 


M 
Mey oor dk 
ss SIND agers ie ee =e oe eer 
0 
+gas [Foe @-m=0, (3) 
M 


where D°. = D..(M); M is the nucleon mass; k = Ez? — M?; o,, are the total 
cross-sections for (pp) and (pp) interactions. 


1) In [9] from an analysis of relation (1) the value /? = 0.082 + 0.008 has been obtained 
what is in good agreement with our result. I am grateful to the author of paper [9] for making 


the result available. 
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With the help of this relation it is possible to estimate the value of the cross- 
sections o, (HZ) in the non-physical region. 

It is extremely important to take into account equalities (2) and (3) in the calcu- 
lations of the values D,.(Z). If one fails to subtract these equalities from the 
corresponding dispersion relations for elastic (xN)- and (NN)-scattering, then 
even a small inaccuracy in the experimental data, because of which the left-hand 
side of equalities (2) and (3) will be different from zero by a very small magnitude 
x, will make a contribution AD, ~«- EH to D..(£). At high energies this contri- 
bution may be much greater than the real values of D. (£). 


2. Elastic Scattering of 7—Mesons on Nucleons 


In Fig. 1 are given the values of D..1) calculated from the dispersion relations. As 
is seen, at high energies the real part of the elastic scattering amplitude is in- 
dependent of the isotopic spin and equals in the lab. system to the constant value 


D,. ~ —0.35- 10-13 [em]. 


Ad 
4 TKS 


Fig. 1. Energy dependence of the real part of the elastic scattering amplitude of a x~-meson on a proton D_ (T) 
(solid curve) and of azt*-meson on a proton D+(T) (dashed line). Lab. system. The values of the kinetic 
energy 7 are givenin the units of Bev; the values of D uy (7) in the units of 107! em. The inaccuracy of the 
values of Di. (7) due to the errors in the values of DY, and f* is AD + = 0.04 


From the expansion of a dispersion relation in powers of E follows the analytical 
expression for this limiting value 


f pp 1 fEdE | 
Di = Zt) + a3 fF m+ tS, w 


i 


") The numerical calculations have been made at the Computing Centre of the USSR Academy 


of Sciences by me f ele ic c iets pa : 
aLhie y means of electronic computer “Strela”, the program being prepared by I. 
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where the constant a is determined b 
See in more details [7]). 


In the c.m.s. Di (EZ) = D°,(E) is equal to 


y the condition o, =o_=o, when EH >a 


z i 
ees LS mh gored Se open [cm] . (5) 


= FE ea ji 


Here hg a Ay = 2.1 - 10-14 em). 


+02 


-02 


a] 3 4 if 
Fig. 2, The same as in Fig. 1, Centre-of-mass system 
30 
LE 
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‘ig. 3. Differential cross section for elastic scattering of m-mesons on nucleons at zero angle in the c.m.8.: 
do/dQ =I*+ D*. The values of the contribution of the real part of the amplitude D? are given. 
Solid curves — (x p); dashed curves — (x p). The values of do/d2 and D? are givenin the units 107 * cm*; 
the values of the kinetic energy T in Bev 
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As is seen, D‘.(H) decreases with the energy increase. The calculated values of 
D. are presented in Fig. 2. In Fig. 3 are given the corresponding values of the 
differential cross section dog/dQ at the angle 0 = 0. 

We se that at HS>1 Bev D%,(Z) <do/dQ and do,/dQ ~ (Eo,/42)? in con- 
trast to the results for (pp)-scattering [5]. 

In the calculations for HZ > 1 Bev the values of the cross-sections from [11] have 
been used. 


3. Elastic Scattering of Nucleons and Antinucleons on Nucleons 


Like for (xN)-interaction, the values of D,. at high energies become constant: 


“9 Le ied ein ee ee eee 1 hea ae 
De = DLT DD) > 39, a di ~~ 0(#) + o-(#)] 


a 
ih ye H o:-a M 
ee — 13; —— 6 
sag | ee i +7Ol+ sas (6) 
M 
where the constant a is determined by the condition o, = o_ =o, when EH >a. 


However, in this case two circumstances made the numeral calculation of more 
difficult: 


a) The unknown values of the cross-section o,(#) in the non-physical region 
Ee [0; 1] enter the dispersion relation. (M is the nucleon mass.) 


b) The experimental values of the (NN)-interaction total cross section is known 
only up to the energies H = 3 Bev. (His the total energy.) 

Nevertheless, in the region of high energies H > M in this case also one can draw 
a conclusion about the magnitude of D, (£). 

The values of o_ experimentally unknown at H > 3 Bev may be obtained by 
extrapolating the values o_, known at lower energies, taking into account at the 
same time that at high energies the equality o_ = o, must be fulfilled. As the 
calculations have shown, the errors which are possible at such an extrapolation 
do not strongly change the magnitude of (pp)-scattering amplitude. 

Taking relation (3) into account, it is easy to see that at H > M with an accuracy 


up to the terms O[x(M/E)], where « < 1, the contribution of the non-physical 
region reduces to the constant: 


M 
ke i} —> const + O[a(M/E)]. (7) 
10) 


The sign and the numeral value of this constant may be determined from the two 
experiments at certain fixed energies H, > M and H,S> M; (e.g. at energies 


8.5 Bev and 6.15 Bev, where the experimental date are in preparation in Dubna 
and Sophy). 
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Fig. 4 represents the results of the numeral calculations of the real part of (pp) 
elastic scattering amplitude 
D,(#) = const + 


1 dE" o,.(E’) 


(ZE" — M?) 


(The values of M; u, EH are given here in Bev; the values of o, in the units of 
10-?’ cm?, the values of D, in the units of 10-3 em). In the calculations for 
H > 1 Bev the values of the cross-sections o4 from [11] have been used. 


(D-C) 


3 
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Fig. 4. Hnergy dependence of the real part of the elastic (pp) scattering amplitude Di(7). Lab. system. The 
values of the kinetic energy 7 are given in Bey; the values of D,(7) — in the units of 1071 cm. Cis a 
constant 


At kinetic energy T =8.5 Bev doq(0)/dQ ~ 130.10-? cm? (in the centre 
of mass system) [5, 6]. 

From Fig. 4 it follows, that at 7’ = 8.5 BeV (D, — C) = 0.1- 10-8 cm. 

Then 


: 2 
CO =0.1-10-8 + - ne (; = = (0.1 + 12.6) + 10-38 om ~ 
a c 


~~ +13-10-% cm. 


Here oc. = 30- 10-7’ cm? [5, 12] 1, =0.99-10-4cm and A= 0.211: 10-4 cm 
are wave lengths in the c. m.s. and in the labor. system. 

We can discriminate between two obtained values of C if we compare these values 
with the experimental magnitudes of do(0)/d2 at 7’ == 8.5 Bev, or if we ana- 
yse the interference of nuclear and Coulomb scattering by very small angles. 
From the point of view of dispersion relation (8)a disagreement between the va- 
ues of the constant C defined at 7’ = 8,5 Bev and 7 =+ 8,5 Bev would imply 
hat the experimental values of do(0)/dQ or o,, used for the calculation of C, are 
1ot accurate. 


4. Phenomenological Analysis of Elastic Scattering of High Energy Particles 


[he amplitude of elastic scattering is represented in this case as an expansion in 
uegendre polynomials 


i xX ee 
A (6) =D (2b + 1) (1 — em) P, (cos 0) (9) 
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with the complex phases 7, which are determined from the comparison with 
experiment. The number of the expansion terms J is determined by the dimension 
R of the spatial region which we intend to investigate. At high energies, where the 
wavelength of the scattered particle is 1<R, N= R/A. 

The complex refraction coefficient x(r) = n(r) + tk(r) is connected with the 
phase 7, by the integral equation 


co 


no) = [ x (Vo? + 8°) ds, (10) 
6 
where 9 =A) 1(l + 3/,) ~ Al (see in more detail [4)]). 
This equation proves to be especially useful in the study of the interaction of high 
energy particles with nuclei, since in this case the analytical form of the function 
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Fig. 5. The cross section for interaction of high energy nucleons with nuclei. Three upper curves are the total 
cross sections o, for the energies 1; 5; 10 Bev (from above to below). Three lower curves are the cross 
sections for inelastic interactions o;, for 1; 5; 10 Bev (from above to below). The values of the cross sec- 
tions are given in the units of 10~*4 em*. A is the atomic weight of nuclei 


x(r) is known from the experiments on quick electron scattering on nuclei. In 
Fig. 5 are given the values of the cross sections for high energy nucleon inter- 
action with nuclei calculated by using relation (10). These values agree well with 
the experimental data known for H = 1.4; 2.8; 4.5 Bev. 

In case of interaction of elementary particles an inverse problem has to be usually 
solved: by the values of 7, known experimentally it is necessary to determine the 
values x(r). In paper [1], to solve this problem a numerical method was applied. 
However, it is possible also to get an analytical solution since equation (10) reduces 
to the integral equation of Abel type by substituting the variable 9? =¢ and 
s? =u —t. The solution of this equation is as follows 1) 


i ee ee 
Qe x dg EN rag (11) 


i re of the numerical calculations of the coefficients are given in papers 


=) Formula (11) was independently obtained also by the workers of the Budapest Physical 
Institute Domoxos and SEBASTYAN [9]. I am grateful to Domokos who communicated the 
results of his investigations unpublished yet. 
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